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accommodate  batches  of  samples  containing  widely  differing 
concentrations  of  test  compound.  A  regression  equation  fup  to  4t.h 
degree  polynomial^  is  calculated  to  describe  the  standard  curve  with 
concentration  as  the  independent  variable  and  absorbance  as  the 
dependent  variable.  Concentrations  in  unknown  samples  are  calculated 
from  absorbance  by  using  a  method  of  successive  approximations.  A 
hard  copy  data  printout  lists  input  data,  concentration  of  test 
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Some  spectrophotometric  procedures  in  analytical  chemistry  do  not 
conform  to  Beer's  law  ( l) .  Examples  of  nonlinearity  between 
concentration  and  absorbance  include  the  thiocyanate  assay  for  iron 
( 2 ),  the  Lowry  method  for  protein  (Si,  and  the  method  for  chloride 
recommended  by  Technicon  (4^.  Nonlinearity  is  especially 
characteristic  of  microbiological  assays  for  nutrients 

Deviation  from  Beer's  law  may  occur  when  the  absorption  band  Cor 
emission  linel  is  not  completely  resolved  Ml.  Most 
spectrophotometers  employ  a  finite  group  of  frequencies.  ’’’he  wider 
the  bandwidth  of  radiation  passed  by  the  filter,  the  greater  the 
apparent  deviation  from  Beer'  law.  The  deviation  is  generally  more 
evident  at  higher  concentrations  and  the  cur/e  bends  toward  the 
concentration  (xl  axis.  Also  discrepancies  are  also  usually  found 
when  the  absorbing  solute  dissociates  or  associates  in  solution, 
because  the  nature  of  the  species  in  solution  will  va^y  with 
concentration.  Scattered  light  from  suspensions,  fluorescence,  or 
stray  light  may  also  cause  deviation  from  Beer's  law  Ml. 

When  an  analyst  arbitrarily  draws  a  straight  line  or  eaHmates  a 
linear  slope  through  standard  curve  points  which  are  not  linear,  but 
curve  toward  the  x  axis,  the  result  is  underestimation  of  test  samples 
with  low  concentrations  and  overestimation  of  the  high  sampler;.  When 
the  nonlinearity  is  severe  enough,  the  analyst  i3  forced  to  ratpr  *n 
the  standard  curve  for  each  test  sample.  This  is  tedious  and  highly 
subject  to  human  error.  The  program  described  in  this  report  was 
created  to  eliminate  the  need  for  visual  reference  t.o  a  standari  curve 
and  to  speed  up  the  processing  of  dots  from  such  analyses.  The 
program  is  also  convenient  for  linear  standard  curves  when  several 
dilutions  per  test  sample  are  assayed. 

PROGRAM  DESIGN 

A  computer  program  wa3  written  to  perform  the  calculations 
necessary  for  analysis  of  curvilinear  spectrophotometric 
absorbance/concentration  data.  The  computer  hardware  and  commercial 
software  necessary  for  implementation  of  the  orogram  are  presented 
under  Materials.  Details  of  the  program  itself,  including  a 
discussion  of  mathematical  algorithms  used,  a  definition  of  derma 
required  for  data  entry,  and  limitations  of  the  program,  are  addressed 
in  Methods.  A  step-by-step  presentation  of  the  program  as  encountered 
by  the  user  during  execution  is  included  under  ’TVER'S  GUTDE. 

Materials 

To  make  the  program  available  to  all  potential  users  al  the 
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'nf!*  i  it-,  wa  s  wri ^en  f'ir  and  implemented  on  the  Institute's 

oen'  ra1  minicomputer,  a  Data  General  Corporation  Rrlipse  G5^0  ( >< 

ft'nrrii  pu rpo:v»  minicomputer  with  51?, 000  words  (16  bits  each'1  of  main 
memory  running  under  control  of  Data  General  Corporation’s  Advanced 
Operating  System  with  1 1  timesharing  support).  The  program  is 
written  in  Fortran  V  (Data  General  Corp.,  Version  6.0?,  a  sunerse4’  of 
Fortran  TV  and  ANSI  Fortran  including  elements  of  TBM  Fortran  TV  and 
Uni  vac  Fortran  V).  For  regression  and  statistical  analysis  of 
standard  curve  data,  the  program  incorporates  a  computational 
subroutine  (RLFOR)  written  in  Fortran  V  by  the  International 
Mathematical  and  Statistical  Libraries  (TMSL'i,  Houston,  Texas  (6). 
This  subroutine  is  part  of  the  statistical  software  package  available 
to  all  users  of  the  Institute’s  central  computing  facility.  Users  can 
gain  access  to  the  program  at  any  of  the  timesharing  terminals  (CRT  or 
printer)  located  throughout  the  Institute.  Hard  copy  output  data  from 
the  program  is  available  to  the  user  on  a  Dataproducts  2250  Line 
Printer  located  at  the  Institute’s  central  computing  facility. 

Methods 

The  following  paragraphs  present  the  mathema+ical  algorithms  used 
in  the  program,  limitations,  and  special  terms  required  for  the  user 
to  understand  and  use  the  program. 

The  mathematical  model  used  to  relate  absorbance  to  concentration 
data  for  purpose  of  generating  a  regression  equation  from  user  input 
standard  curve  data  is  the  simple  polynomial 

y  =  a  +  .a  X1  t  .  +a  x° 

o  1  n 

where,  x  =  standard  concentration  (independent  variable) 
y  -  observed  absorbance  (dependent  variable) 
n  =  maximum  degree  of  polynomial  (discussed  below) 

This  model  allows  for  the  calculation  of  a  regression  equation  which 
can  vary  from  a  linear  model  (where  n  =  1)  up  to  a  poTynomia1  of 

maximum  degree  n.  The  maximum  degree  of  the  regression  equation 
allowed  by  the  program  has  been  arbitrarily  limited  to  four.  We  have 
found  that  polynomial  equations  up  to  a  maximum  degree  of  four  (a 
quart’ c  equation)  can  he  used  to  model  even  those  assays  which  exhibit 
extreme  nonlinearity  (e.g.,  microbiologial  nutrient  assays). 

The  TMRL  subroutine  (6)  included  in  the  program  performs  ’■he 
computations  required  for  a  standard  univariate  curvilinear  regression 
analysis.  Orthogonal  polynomals  are  used  by  the  subroutine  to  compute 
the  regression  polynomial  (7,R).  Tn  addition  to  calculating  the 
regression  coefficients  of  the  polynomial,  the  subroutine  also 
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provides  the  usual  analysis  of  variance  table  information  along  with 
other  statistical  inferences  for  each  model  parameter.  The  fitted 
model  returned  by  the  subroutine  is 


where : 


yi 


yi 


-  a  +  a4x.  +  *  ’  *  * x?  +  p. 

o  1  i  2  i  pi  i 

till 

=  predicted  absorbance  for  i  '  data  point 
(dependent  variable! 


x. 

i 


standard  concentration  for  i  data  point 
(independent  variable! 


p  =  maximum  decree  of  polynomial 

(limited  to  four  in  this  program! 


a^  =  computed  regression  coefficients 

i.  =  1  , ? . .  N  (number  of  data  points! 

e^  =  population  random  error 

The  population  random  errors  corresponding  to  e.  a"®  assumed  to  be 
uncorrelated,  to  have  means  0*"  zero,  and  have  variances  6~  '  .  billing 
execution  of  the  program  the  user  has  the  choice  of  computing  either  a 
polynomial  regression  model  or  a  linear  regression  model  (;n  cases 
where  it  is  known  that  the  relationship  between  absorbance  and 
concentration  is  linear  under  given  assay  conditions.! 


Because  of  the  possibility  that  the  generated  polynomial  will  no* 
be  monotonic  over  the  input  concentration  range,  a  subroutine  has  been 
incorporated  into  the  program  to  detect  this  condition.  Tf  the 
generated  polynomial  is  not  monotonic  over  the  input  concentration 
range,  the  program  will  abort  and  the  user  will  be  instructed  *0  che^k 
data  for  accuracy  or  prepare  a  new  standard  curve  if  necessary. 


In  order  to  facilitate  the  computation  of  sample  concent ra *' on 
data  from  the  calculated  regression  model  and  also  to  a1 low  the 
analyst  a  certain  degree  of  latitude  in  preparing  standard  curves,  the 
concepts  of  "standard  volume"  and  "standard  equivalen*  volume"  have 
been  used  (see  Appendix  A,  Table  1  for  definitions!.  Tn  preparing  *he 
standard  curve  for  a  given  assay,  it  will  be  assumed  that  the  analyst 
has  prepared  one  solution  of  known  concentration  (representing  the 
maximum  concentration  in  -the  curve!  and  from  this  solution  will 
prepare  the  standard  curve.  For  purposes  of  i l1 nstr«tf on,  a 
hypothetical  example  is  given  and  the  following  assay  parameters 
assumed : 
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a.  Total  volume  of  standard  (including  diluent."'  =  2.0  ml 

b.  Number  of  points  in  standard  curve  =  11. 

The  analyst  woul  i  prepare  his  standard  cirve  in  the  foil' wing  manner. 


STANDARD 

VOLUME  OF 

SS  TOTAL 

STANDARD  EQUIVALENT 

VOLUME* 

DILUENT 

VOLUME 

VOLUME 

C ml  1 

Cml  1 

(mil 

Cml  1 

Co.o) 

2.0 

2.0 

0.0 

Co.2") 

1  .fi 

2.0 

0. 1 

(0. 4 1 

1.6 

2.0 

0.2 

Co.6l 

1  .4 

2.0 

0.1 

Co.nl 

1.2 

2.0 

0.4 

C 1  .O'! 

1  .0 

2.0 

0.6 

(1  .2") 

O.R 

2.0 

o.f, 

(1 .4"! 

0.6 

T.O 

0."' 

Cl  .hi 

0.4 

2.0 

O.R 

Cl.nl 

0.2 

2.0 

O.R 

C?.0l 

0.0 

2.0 

1  .0 

The  user  enters  only  the  "Standard  Volume"  ("note!  by  an  ssterisV 
shovel  along  with  its  correspond ing  absorbance  during  execution 
of  the  program.  The  "standard  volume"  is  scaled  by  the  program 
over  the  range:  0.0  <_  STANDARD  EOUTVALENT  VOLUME  1.0  for 
computational  purposes  within  the  program  and  for  graphicpi 
presentation  of  the  generated  regression  curve  in  the  hardcopy 
data  output  of  the  program.  Because  of  this  scaling  feature,  *-he 
user  is  free  to  choose  the  assay  volume  to  suit  the  test  and/or 
instrumentation  requi remenis.  The  maximum  absorbance  accepted  by 
the  program  in  its  current  form  is  2.0.  Tf  transmittance  values 
a me  determined  instead  of  absorbance  units,  they  must  be 
•’onverted  by  the  user  into  absorbance  uni  t.s  before  the  data  can 
be  entered  in  the  program.  A  total  of  MO  standard  curve  data 
points  ‘  i ,  e  •  ,  absorbance,  standard  volume  pairs’!  can  be  input. 

Upon  successful  completion  of  the  regression  equation  ' see 
Appendix  A,  Table  2  for  erm o1'  code  descriptions!,  the  program 
accepts  certain  alphameric  information  used  in  the  hardcopy  data 
output.  Thi3  information  consists  of  l)  name  of  assay,  2)  type 

of  sample,  Ml  date  of  assay,  and  Al  name  of  investigato”. 

Complete  instructions  for  entering  this  information  are  provided 
under  USER'S  OUTRE.  Standard  concentration  and  sample  dilution 
factor  information  is  also  entered  at  this  tim®,  along  with  +h<» 

number  of  dilutions  par  sample  and  the  volume  of  the  sample  at 

each  dilution.  The  program  win  accept  up  to  A  dilutions  per 
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sample.  Finally,  fhe  sample  number  and  the  corresponding 
measured  absorbance(s)  (if  more  than  one  sample  volume  was 
prepared)  for  each  sample  are  entered. 

The  mathematical  technique  used  in  this  program  to  compute 
concentration  from  measured  absorbance  is  a  numerical  iterative 
procedure  commonly  referred  to  as  the  method  of  successive 
approximations.  Because  the  standard  curve  regression  equation 
was  generated  with  absorbance  (dependent  variable)  as  a  function 
of  concentration  (independent  variable),  the  equation  must  be 
either  solved  analytically  for  concentration  as  a  function  of 
absorbance  or  numerical  techniques  must  be  applied  to  arrive  at  a 
solution.  Although  analytical  methods  are  available  for  the 
closed  solution  of  equations  up  to  the  third  degree,  no  generally 
applicable  technique  exists  for  the  solution  of  quart.ie  and 
higher  degree  equations  (0,10).  Iterative  techniques  approximate 
the  root(s)  of  equations  to  some  prescribed  degree  of  accuracy. 
The  subroutine  employed  in  this  ^program  computes  the  unknown 
concentration  to  +  1.0  x  10”  standard  equivalent  volume.  A 

complete  description  of  this  technique  is  presented  in  Appendix 
B,  under  Subroutine  'BSECT'. 

After  the  unknown  "standard  equivalent  vo1ume"  has  b»en 

determined  from  the  measured  absorbance,  +he  following 

relationship  is  used  to  compute 

the  actual  concentration  of  subs+anee  in  the  samplers'*. 

SUBSTANCE]  =  (STD.BQ.VOLVSS.TOT.VOLVSTD.CONC.  Vi  /SU.VOLVnF) 

where : 

SUBSTANCE  =  concentration  (or  amount/unit.  weight)  of  test 

substance  in  original  sample. 

STB. EQ. VOL.  =  Standard  equivalent,  volume  ^ calculated  by  BSECT 

subroutine) . 

SS. TOT. VOL.  =  sample  and  standard  total  volume. 

STD.CONC.  =  concentration  of  standard  used  to  prepare  curve. 
SD.VOL.  =  sample  dilution  volume. 

DF  =  dilution  factor. 

(See  Appendix  A,  Table  1  for  further  explanations  of  +  he 
abbreviations . ) 


The  preceding  calculation  is  based  on  RS  total  volume  and  is  therefore 
independent,  of  total  reaction  volume.  ?oT*  purposes  of  this  program, 
the  amounts  of  any  additional  reagents  do  not  affect  the  col rul n*i ons . 

The  program  also  provides  the  option  of  calculating  more  than  one 
set  of  sample  data  (with  different  dilution  factors,  sample  volumes, 
ete.l  from  the  same  standard  curve  regression  equation.  A  complete 
listing  of  the  program  is  provided  in  Appendix  B. 


USER'S  GUTDE 


General 

Before  giving  a  detailed  description  of  the  program  from  the 
user's  point  of  view,  two  points  must  he  discussed.  The  first  deals 
with  limitations  of  the  program  and  the  second  iB  related  to 
provisions  made  for  'unknown'  sample  data  entry. 

This  program  provides  the  analyst  with  a  simple  and  accurate 
computational  method  for  evaluating  spectrophotometric  assay  data  as 
an  alternative  to  graphical  techniques  which  sre  usually  tedious  and 
error-prone.  The  program  generates  a  linear  or  polynomial  regression 
equation  from  input  standard  curve  data  points.  This  equation  is  then 
used  as  a  statistical  model  to  estimate  unknown  sample  concentration 
data  from  known  absorbances.  The  program  should  not  be  used  for 
purposes  of  inferring  information  about  the  nature  of  the  chemical 
reaction  involved  in  a  particular  chemical  or  microbiol ogical  assay 
(e.g.  reaction  mechanisms  or  kinetics). 

Provision  has  been  made  in  the  program  to  allow  entry  of  un  to 
four  dilutions  per  'unknown'  test  sample.  There  are  several  reasons 
for  including  this  feature.  First,  interfering  substances  can  often 
be  detected  by  assaying  the  test  sample  at  several  differen+ 
dilutions.  This  problem  is  common  in  microbiological  assays  and  is 
referred  to  as  "upward"  or  "downward"  drift,  depending  on  whether  the 
nonspecific  interference  is  stimulatory  or  inhibitory.  Thus  for  some 
assays,  this  feature  could  be  considered  as  a  monitor  of  the  validity 
of  the  results  for  each  sample.  Secondly,  with  several  dilutions,  the 
analyst  can  handle  samples  with  widely  different  concentrations  of 
test  substance  without  being  forced  to  repeat  samples  which  fall  out 
of  the  range  of  the  standard  curve.  It  is  often  time-effective 
( particularly  if  automated  pipetting  equipment  is  available)  to 
include  extra  dilutions  and  to  disregard  outlying  values.  Finally, 
the  feature  increases  the  statistical  power  of  the  resultant,  data  by 
increasing  the  number  of  replicates. 

Specific  User  Instructions 

In  the  complete  operational  description  of  fhe  program, 
capitalized  sequences  denote  program  output  displayed  on  the  user's 
terminal,  underlined  sequences  denote  prompts  requiring  user  date 
input,  and  the  appropriate  user  response  is  enclosed  in  parentheses. 
The  program  is  outlined  in  the  following  nine  steps. 


STEP  1.^  The  user  accesses  the  program  at  a  timesharing  terminal  by 
entering  the  CLT  command  'EXECUTE  ASSAY'.  A  brief  introduction  to  the 
program  will  be  displayed  on  the  terminal.  The  user  presses  <RETURN> 
to  resume  program  execution. 


STEP  2.)  Instructions  are  given  for  entering  standard  curve  data 
points,  then  the  following  series  of  prompts  are  displayed  for  data 
point  entry. 


DATA  POINT  NO.:  'n' 


ENTER  VOLUME  OF  STANDARD  ( ' IND  VARIABLE'): 

~  TUser  enters  volume  of  standard  used  fo~ 
given  data  point.) 

ENTER  MEASURED  ABSORBANCE  FOR  ABOVE  VOLUME  ('PEP'  VARIABLE): 

(User  enters  absorbance.  Allowed 
absorbance  range  is  0.0  to  2.0.) 

EXAMINE  ABOVE  ENTRIES.  IF  YOU  WISH  TO  CORRECT 

EITHER  OR  BOTH  VALUE(S),  ENTER  T  BELOW. 

. ENTER  T  TO  CORRECT  DATA  POINT  ENTRY. 

. ENTER  '2'  TO  CONTINUE  DATA  POINT  ENTRY. 

. ENTER  '999'  TO  TERMINATE  DATA  POINT  ENTRY. 

ENTER  CHOICE:  (User  enters  '1'  to  correct  data  point 

entry  just  made,  to  enter  another 

data  point,  or  '999'  to  terminate  data 
point  entry.) 


At  least  five  data  points  are  recommended  for  linear  standard  curves 
and  at  least  eight  data  points  are  recommended  for  nonlinear  standard 
curves.  A  maximum  of  thirty  data  points  may  be  entered. 
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STEP  ■?.)  After  all  standard  curve  data  points  have  been  entered,  the 
program  outputs  to  the  terminal  the  following  message. 


ENTER  TYPE  OP  REGRESSION  CURVE  DESIRED  —  LINEAR  OR  POLYNOMIAL 


ENTER  *1‘  FOR  LINEAR  REGRESSION  CURVE. 
ENTER  '2'  FOR  POLYNOMIAL  REGRESSION  CURVE. 


ENTER  CHOICE:  (User  enters  '1'  if  the  assay  is  known  to 

be  linear,  or  '2'  if  a  polynomial 
regression  curve  is  desired.^ 


STEP  4.)  Program  displays  the  following  message  at  the  terminal: 
CALCULATING  POLYNOMIAL  REGRESSION  EQUATION. 


If  any  error  conditions  are  detected  by  the  regression  subroutine,  a 
message  indicating  the  fact  will  appear  at  the  terminal.  The  most 
likely  cause  for  an  error  to  be  detected  is  an  insufficient  number  of 
data  points  for  the  requested  type  of  regression  equation.  If  an 
error  occurs,  the  program  must  be  aborted  and  the  cause  corrected 
before  the  user  can  continue.  Possible  error  conditions  are  listed  in 
Appendix  A,  Table  2. 


If  the  regression  equation  was  successfully  computed,  the  following 
message  appears  at  the  terminal: 


POLYMONIAL  REGRESSION  EQUATION  OF  DEGREE  V  HAS  BEEN  CALCULATED. 

STEP  5.)  The  user  is  next  requested  to  enter  information  to  annotate 
output  data.  The  following  prompts  are  displayed: 


ENTER  TYPE  OF  ASSAY: 

User  enters  discription  of  assay. 
(Maximum,  ^2  alphameric  characters.'' 


9 


ENTER  DATE  OF  ASSAY: 

(User  enters  the  date  assay  was 
conducted.  Maximum  of  12  alphameric 
characters  may  he  entered.'' 

ENTER  TYPE  OF  SAMPLE: 

(User  enters  description  of  sample. 
Maximum  of  12  alphameric  characters  may 
be  entered.’' 


ENTER  NAME  OF  INVESTIGATOR; 

(User  enters  name  of  person  conducting 
assay.  Maximum  of  12  alphameric 
characters  may  be  entered.'' 

STEP  Information  required  to  calculate  'unknowns’  is  now  entered. 
These  data  consist  of  dilution  factor  information,  concentration  of 
standard,  the  number  of  sample  dilutions  and  the  volume  of  each  sample 
d  ? 1 ut ion. 


STANDARD  CONCENTRAT CON  DATA: 


ENTER  NUMERICAL  PORTION  OF  STANDARD  CONCENTRATION: 

(User  enters  numerical  portion  of 
'standard'  concentration."' 


. ENTER  UNITS  OF  STANDARD: 

(User  enters  units  of  standard  (e.g., 
ng/ml,  mg/ml,  etc.  A  maximum  of  6 
characters  may  be  entered.) 

ENTER  INITIAL  DILUTION  F.tCTOR  OF  SAMPLE  MATERIAL: 

(User  enters  initial  dilution  factor  of 
sample  material.) 
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Tf  initial  dilution  factor  is  greater  than  1000,  the  user  will  be 
requested  to  re-enter  the  units  of  the  standard.  For  example,  if  the 
original  units  were  'ng/ml'  then  ' ug/ml'  should  be  entered,  or  if  the 
origins''  units  were  'mg/l‘  then  g/l  should  be  entered. 


ENTER  NUMBER  OF  SAMPLE  DILUTIONS  PER  SAMPLE  (MTN.  1^  MAX.  4): 

(User  enters  number  of  sample  dilutions 
made  for  esch  sample.  If  only  one,  a  '1* 
must  be  entered.  Maximum  that  can  be 
entered  is  ' 4 ' . 1 


ENTER  'MLS  OF  SAMPLE’  AT  DILUTION  'n' ; 

(User  enters  the  number  of  milliliters  of 
sample  used  at  each  dilution  volume.'' 


STEP  s.l  Unknown  absorbances  snd  sample  numbers  are  now  entered  in+o 
the  program.  A  brief  description  of  the  sample  entry  procedure  is 
displayed  at  the  terminal,  then  the  following  prompts  are  given. 


ENTER  SAMPLE  NUMBER: 


(User  enters  a  numerical  sample  number. 
The  range  the  sample  number  can  assume  is 
from  i  to  09999.) 

ENTER  ABSORBANCE  AT  SAMPLE  PTLUTION  ’n‘ : 

(User  enters  measured  absorbance  at. 
indicated  sample  dilution.) 

EXAMINE  ABOVE  DATA.  IF  YOU  WISH  TO  CORRECT  ANY 
ENTRY,  ALL  DATA  FOR  SAMPLE  (n)  MUST  BE  RE-ENTERED. 

. ENTER  T  TO  CORRECT  DATA. 

. ENTER  ’2'  TO  CONTINUE  WITH  DATA  ENTRY. 

. ENTER  '999'  TO  TERMINATE  DATA  ENTRY. 
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ENTER  CHOICE: 


(User  enters  *  1'  i  f  he  wishes  to  eorreo4' 
a  date  entry  ,iust  made,  to  erter  data 

for  another  sample,  or  '009*  to  terminate 
data  entry. 1 


The  user  may  enter  absorbance  dsts  for  up  to  fifty  (50)  samples. 


STEP  p.)  The  user  is  now  given  the  option  of  entering  another  set  of 
data  using  the  same  regression  our/e.  For  example,  if  the  user  had 
dsts  for  more  than  r>0  samples,  he  could  use  this  feature  to  avoid  the 
necessity  of  re-calculating  the  standard  curve.  This  feature  is  also 
helpful  if  some  samples  had  different  initial  treatment  (e.g., 
different  dilution  factor). 


HO  YOU  WISH  TO  ENTER  MORE  DATA  USING  SAME  STANDARD  CURVED 


. ENTER  ' 1 '  TO  COMPUTE  ANOTHER  DATA  SET . 

. ENTER  '?'  TO  EXIT  PROGRAM. 

ENTER  CHOICE: 

Miser  enters  '1'  if  he  wishes  to  compute  another  set  of  data  using 

same  regression  curve,  or  '2' to  exit 
program. ) 


STEP  9.’'  At.  this  point  the  following  message  is  output  to  users 
•■arnMnal  and  program  execution  ceases. 


OUTPUT  DATA  FThE  'ASSAY. DATA'  CREATED.  END  OF  PROGRAM. 


An  output  data  file  named  'ASSAY. DATA'  was  created  by  the  program 
in  the  user's  directory.  A  hardcopy  printout  of  this  file  can  be 
obtained  at  the  line  printer  located  in  the  Institute's  central 
computing  facility  by  entering  the  CLI  command  'QPRTNT  ASSAY. DATA' . 

Hardcopy  output  provided  by  the  program  to  the  user  consists  of 
regression  equation  statistical  data  (analysis  of  variance  tablp 
information  along  with  statistical  inferences  pertaining  to  computed 


regression  coefficients),  input  standard  curve  dats  points,  and 
predicted  values  of  absorbance  ( dependent  variable)  calculated  using 
the  regression  model.  A  low  resolution  graph  of  both  the  input 
standard  curve  data  points  and  the  computed  regression  equation  is 
provided  for  visual  inspection  of  the  curve.  Sample  dat8  output 
includes  assay  parameters,  input  sample  absorbances  along  with 
corresponding  actual  concentration  of  test  substance  in  the  sample 
material.  If  more  than  one  dilution  per  test  sample  was  performed, 
the  mean  concentration  is  calculated.  The  output  includes  not.  only 
the  mean  of  all  calculated  values  for  esch  sample,  but  also  an 
‘adjusted  mean'.  The  choice  of  the  range  (mean  _+  10^)  for  accepting 
values  for  the  adjusted  mean  was  purely  arbitrary.  The  purpose  of 
this  feature  is  not  only  to  eliminate  spurious  values,  but  also  to 
allow  the  analyst  to  ascertain  the  variability  of  the  assay.  By 
scanning  one  column  in  the  output,  one  can  quickly  determine  how  many 
(if  any)  dilutions  yielded  values  which  were  more  than  10^  above  or 
below  the  mean.  Examples  of  data  output  for  both  a  linear  and 
nonlinear  assay  are  given  in  Appendix  C. 

CONCLUSIONS 

The  computer  program  described  in  this  report  provides  a  simple 
and  accurate  method  for  evaluating  spectrophotometric  assay  data.  It 
is  a  general  purpose  program  which  will  process  any  analytical  data 
that  must  be  referred  to  a  standard  curve.  It  will  accommodate  up  to 
4  dilutions  per  sample  and  thus  can  handle  batches  of  analyses  with 
widely  differing  concentrations  of  test  substance.  The  program  is 
available  to  anyone  at  this  Institute. 

RECOMMENDATIONS 


None . 
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TABLE  1  -  DEFINITION  OF  TERMS  USED  IN  PROGRAM 


STANDARD  VOLUME:  Actual  volume  (in  milliliters)  of  standard  used  to 
prepare  each  individual  point  in  standard  curve. 

SS.  TOTAL  VOLUME:  Standard  or  sample  total  volume  (volume  of  standard 
or  sample  plus  diluent).  This  vslue  must  equal  maximum  volume  of 
standard. 


STANDARD  EQUIVALENT  VOLUME:  Input  standard  volumes  scaled  over  the 
range  of  0.0  to  1.0  inclusive.  The  standard  volumes  are  scaled 
over  this  range  for  computational  purposes  and  graphical 
presentation  of  the  standard  curve  in  the  program  output  data. 


DILUTION  FACTOR:  Initial  dilution  factor  of  sample  material  which 
occurs  during  sample  preparation  including  dilution  necessary  to 
bring  absorbance  of  test  substance  into  range  of  the  standard 
curve. 


STANDARD  CONCENTRATION :  Concentration  of  test  substance  in  the 
standard  used  to  prepare  the  standard  curve. 


SAMPLE  DILUTION  VOLUME:  Actual  volume  of  sample  (in  milliliters) 
included  in  reaction  mixture.  Thi3  volume  plus  diluent  must  equal 
the  SS. TOTAL  VOLUME. 


MEAN  CONCENTRATION :  Computed  average  concentration  of  test  substance 
in  sample  if  more  than  one  sample  dilution  volume  analyzed. 


ADJUSTED  MEAN  CONCENTRATION:  Computed  average  concentration  of  test 
substance  in  sample  excluding  those  values  exceeding  ♦_  of  the 
computed  mean  concentration  if  more  than  one  sample  dilution 
volume  anelyzed. 
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TABLE  2  -  PROGRAM  ERROR  CODES 


STANDARD  CURVE  ENTRY :  All  errors  detected  by  the  program  which  are 

related  to  standard  curve  data  point  entry  are  fatal  to  program 
execution-  The  user  must  correct  the  problem  in  the  data  and 
then  access  the  program  again.  Errors  at  this  point  in  the 
program  are  always  caused  by  failure  to  enter  at  least  three  data 
points  for  a  linear  regression  model,  or  a  minimum  of  five  data 
points  for  a  polynomial  regression  model.  If  an  error  is 
detected,  the  following  message  will  be  displayed  on  the  users 
terminal:  STOP  "INADEQUATE  NUMBER  OF  DATA  POINTS". 


SAMPLE  HARDCOPY  DATA  OUTPUT ;  The  following  error  conditions  may  be 
detected  by  the  program  during  computation  of  output  data. 

1. )  If  measured  absorbance  of  test  substance  is  outside  the  range 
of  the  computed  standard  curve  regression  model,  an  error  code  of 
'-111'  will  appear  in  the  output  data  of  the  program  in  the 
appropriate  column.  To  correct  this  error  the  user  must  dilute  sample 
material  to  bring  concentration  of  the  test  substance  into  range  of 
the  standard  curve. 

2. )  if  the  concentration  of  te3t  substance  at  all  dilution 
volumes  for  a  given  sample  exceeded  +_  10^  of  the  computed  mean 
concentration,  an  error  code  of  '-NNV  will  appear  in  the  'ADJUSTED 
MEAN  CONCENTRATION’  column  of  the  output  data. 
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Q  ********************************************************************** 

c.  *  * 

C  *  ASSAY  DATA  ANALYSIS  PROGRAM  -  PROGRAM  NAME:  "ASSAY"  * 

C  *  * 

C  *  WRITTEN  -  MAY  1981  -  J..7.  KNUDSEN,  LAIR,  DTV.OF  RESEARCH  * 

C  *  SUPPORT,  ANALYTICAL  CHEMISTRY  CROUP  * 

C  *  * 

C  *  LANGUAGE  -  DATA  GENERAL  FORTRAN  V  * 

C  *  * 

C  *  PURPOSE  * 

C  * -  * 

C  *  PROGRAM  TO  PROCESS  SPECTOOPHOTOMETRIC  ANALYTICAL  DATA  ASSOCIATED  * 

C  *  WITH  CURVILINEAR  ABSORBANCE/CONCENTRATION  RELATIONSHIPS  * 

C  *  * 

C  *  PROGRAM  SYMBOL  DEFINITION  TABLE  (EXCLUDING  LOOP  CONTROL  INDEXES)  * 

C  * -  * 

C  *  INTEGER  VARIABLES  * 

C  * -  * 

C  *  IX  -  ARGUMENT  SUBRT  'RLFOR' .INPUT  ROW  DIM  OF  MATRIX  'XYW'  * 

C  *  N  -  ARGUMENT  SUBRT  ' RLFOR INPUT  NUMBER  OF  DATA  POINTS  (STD  CURVE)  * 

C  *  MDP  -  ARGUMENT  SUBRT  'RLFOR ' .CONTROL  VECTOR  OF  LENGTH  ^  * 

C  *  IB  -  ARGUMENT  SUBRT  'RLFOR' .INPUT  ROW  DIMENSION  OF  MATRIX  'B*  * 

C  *  IP  -  ARGUMENT  SUBRT  'RLFOR' .INPUT  ROW  DIMENSION  OF  MATRIX  'PRED'  * 

C  *  TER  -  ARGUMENT  SUBRT  'RLFOR' .ERROR  PARAMETER  * 

C  *  I ERROR  -  ARGUMENT  SUBRT  'MONOT' .ERROR  PARAMETER  * 

C  *  I ERR  -  ARGUMENT  SUBRT  'BSECT' .ERROR  PARAMETER  * 

C  *  NVOL  -  MAX  NUMBER  OF  DILUTION  VOLUMES  FOR  GIVEN  ASSAY  * 

C  *  SNUM  -  INPUT  VARIABLE.  TEMPORARILY  STORE  SAMPLE  NUMBER  * 

C  *  SNUMBER  -  VECTOR  OF  MAX  LENGTH  SO  TO  HOLD  SAMPLE  NUMBFRS  * 

C  *  GRAPH  -  TWO  DIMENSIONAL  MATRIX  TO  HOLD  LOW  RESOLUTION  GRAPH  * 

C  *  LABEL!  -  VECTOR  OF  LENGTH  51  TO  HOLD  VERT  AXIS  LAP EL  OF  'GRAPH'  * 

C  *  A  -  TEST  INTEGER  AND  ARGUMENT  OF  SUBRT  'PLOT?'  * 

C  *  MTEST  -  TEST  INTEGER . DETERMINES  IF  MORE  THAN  ONE  SAMPLE  DATA  * 

C  *  ANALYSIS  WILL  BE  PERFORMED  USING  SAME  'STANDARD  CURVE'  * 

C  *  ITEST  -  TEST  INTEGER. CONTROLS  TYPE  OF  REGRESSION  CURVE  CALCULATED  * 

C  *  * 

C  *  REAL  VARIABLES  * 

C  * * 

C  *  ABSORBANCE  -  TWO  DIM  MATRIX  TO  HOLD  INPUT  SAMPLE  ABSORBANCES  * 

C  *  CONC  -  WO  DIM  MATRIX  TO  HOLD  DATA  FOR  SAMPLE  'UNKNOWN'  CONC  * 

C  *  LABEL?  -  VECTOR  OF  LENGTH  51  TO  HOLD  VERT  AXIS  LABEL  OF  'GRAPH'  * 

C  *  MCPLUS  -  VARIABLE  TO  HOLD  'MEAN  CONC  +  (MEAN  CONC) *o  DATA  * 

C  *  MCMINUS  -  VARIABLE  TO  HOLD  'MEAN  CONC  +  (MEAN  CONC) *Q  DATA  * 

C  *  STAND  -  VARIABLE  TO  HOLD  NUMERICAL  PORTION  OF  STANDARD  CONC  * 

C  *  DILFAC  -  VARIABLE  TO  HOLD  SAMPLE  DILUTION  FACTOR  * 

C  *  SVOL  -  VARIABLE  TO  TEMPORARILY  HOLD  INPUT  SAMPLE  VOLUMES  * 

C  *  SVOLUME  -  VECTOR  OF  LENGTH  4  TO  HOLD  INPUT  SAMPLE  VOLUMES  * 

C  *  DAT  -  WO  DIM  MATRIX  TO  HOLD  STANDARD  CURVE  INPUT  DATA  * 

C  *  MAXVOL  -  VARIABLE  TO  HOLD  MAXIMUM  STANDARD  CURVE  INPUT  VOLUME  * 

C  *  MAXABS  -  VARIABLE  TO  HOLD  MAXIMUM  STANDARD  CURVE  INPUT  ABSORBANCE  * 
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r  *  VALUE  -  VARIABLE  TO  TEMPORARILY  HOLD  INFOT  STD  CURVE  TNPUT  DATA 

r  *  DTV  -  VARIABLE  TO  DIVIDE  SAMPLE  CON C  BY  1000  IF  1 DT LFAC '  >  1000. 

C  *  Q  -  VARIABLE  SET  TO  . 1  TO  ARBITRARILY  CONTROL  EXCLUSION  OF  SAMPLE 
C  *  CONCENTRATION  FROM  MEAN  ATXTJSTED  CONCENTRATION  TABULATION 

C  * 

C  *  DOUBLE  PRECISION  VARIABLES 
C  * - 

C  *  XYW  -  ARGUMENT  SURRT  'RLFOR'.  WO  DIM  MATRIX  (DATA  POINT  I/O) 

C  *  RSO  -  ARGUMENT  SUBRT  'RLFOR' .CONTROL  DEGREE  OF  FITTED  MODEL 

C  *  ALBP  -  ARGUMENT  SUBRT  'RLFOR '  .VECTOR  OF  LENGTH  2  CONTAINING  RISK 
C  *  LEVELS 

C  *  ANOVA  -  ARGUMENT  SUBRT  'RLFOR' .OUTPUT  VECTOR  OF  LENGTH  13 

C  *  B  -  ARGUMENT  SUBRT  'RLFOR' .WO  DIM  MATRIX  (OUTPUT  OF  REGRESSION 

C  *  EQUATION  COEFFTCIEOTS  AND  OTHER  STATISTICAL  DATA 

C  *  WK  -  ARGUMENT  SUBRT  'RLFOR '.WORK  VECTOR 

C  *  PRED  -  ARGUMENT  SUBRT  'RLFOR' .OUTPUT  'N*  BY  6  MATRIX 

C  *  XTOL  -  ARGUMENT  SUBRT  ' BSECT' .  SET  EQUAL  TO  TOLERANCE  REOUIRED 

C  *  BY  'BSECT'  SUBROUTINE 

C  *  TEMPI  -  TEMPORARY  STORAGE  VARIABLE 

C  *  TEMP?  -  TEMPORARY  STORAGE  VARIABLE 

C  *  XLO  -  ARGUMENT  SUBRTS  'BSECT' , 'MONOT' , 'PLOT?' .LOWER  RANGE  OF 
C  *  'STANDARD  FOUI VALENT  VOLUME’ 

C  *  XHI  -  ARGUMENT  SUBRTS  'BSECT' , 'MONOT' ,' PLOT? ' .UPPER  RANGE  OF 
C  *  'STANDARD  EOU I VALENT  VOLUME' 

C  *  HD1  -  VECTOR  OF  LENGTH  8  CONTAINING  ALPHAMERIC  ASSAY  DATA 

C  *  HD?  -  VECTOR  OF  LENGTH  6  CONTAINING  ALPHAMERIC  ASSAY  DATA 

C  *  HD 3  -  VECTOR  OF  LENGTH  8  CONTAINING  ALPHAMERIC  ASSAY  DATA 

C  *  HD4  -  VECTOR  OF  LENGTH  8  CONTAINING  ALPHAMERIC  ASSAY  DATA 

C  *  HD5  -  VARIABLE  CONTAINING  ALPHAMERIC  ASSAY  DATA 
C  *  COEFF  -  VECTOR  OF  LENGTH  6  CONTAINING  REGRESSION  COEFFICIENTS 
C  *  SARS  -  ARGUMENT  SUBRT  ' BSECT '. TRANSFER  SAMPLE  ABSORBANCE  TO 
C  *  SUBRT  'BSECT' 


★ 

* 

* 

* 

* 

* 

★ 

* 

* 

★ 

★ 

* 

* 

* 

★ 

* 

* 

★ 

★ 

* 

* 

★ 

* 

★ 

* 

★ 

★ 

ir 

it 

* 

* 

•k 


c  *  * 

Q  ********************************************************************** 


C  PROGRAM  DATA  TYPE  STATEMENTS 

INTEGER  IX,N,MDP(3) , IB, IP, IER, TERROR, NVOL,SNUM,IERR, 

1  SNUMBER(SO) , GRAPH (0:50,0: 100) ,LABEL1 (51) ,A, 

1  ITEST , MTEST , NTEST 

REAL  ABSORBANCE (50, 4) ,CONC(50,8) ,LABEL2(51) ,MCPLUS,MCMINUS, 

1  STAND, DILFAC , SVOL, SVOLUME (4 ) ,DAT(30,3) ,MAXVOL,MAXABS, 

l  VALUE, DIV,Q 

DOUBLE  PRECISION  XYW(30,7) ,RSQ, ALBP(2) ,ANOVA (13) ,B(7 , 12) ,WK (100) , 
l  PRED(30,fi) , XTOL, TEMPI, TEMP2, XLO, XHI, HD8, 

1  COEFF (8) ,HD1 (4 ) ,HD2(4) ,HD3(4) ,HD4(4) ,HD5,SABS,X 
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C  PROGRAM  PARAMETERS 


PARAMETER  XHI=1.0,  XLO=1.0F.-R,  XTOL=].OE-S,  0=1  .OE-1 
C  PROGRAM  SYMBOL  INITIALIZATION'  STATEMENTS 


DATA 

(LABELl(I) ,1=1 

,51)/' 

*«••«!«« 
9  9  9  9 

till 
t  9 

1  \  «  1  1  •  t 

9  9  9  9 

1  V  * 

t  1  1  q  1  «  rp  I  f 

9  9  9  1  9 

A' ,'N' 

9 

O'  /A  VR  VD"  / 

•  »V«  1 

9  V  / 

O' ,’L' ,'U' ,'M' 

1,'EV 

VEVOVU' 

#'IV 

V' 

,'A  VL1  ,*E  VN1 

Iml  •  • 

t  1  9 

1  «  1  •  •  1  • 
9  9  9  9 

l  V  • 

•  1  1  1  1  1  t 

t  f  r  ¥ 

1  1  1 
9 

9 

1  / ,  (LABFL2  (.7),  ,7 

=1,51)/' 

n.ov  '  , '  ' 

1/  V 

»  *  0  1  1  1  «  1 
/  V  •  J  9  9 

f  1  « 

9 

9 

i  » n  o »  »  i  »  * 

9  1  •  *  9  9 

111! 
9  9 

l,'  V 

*  •  n  A  •  *  •  • 

9  ‘  '  •  4  9  9 

1  1  » 
9 

9 

i  in  ci  i  •  i  l 
9  3  9  9 

1  1  1  1 
9  9 

•  n  fi*  '  »  •  • 

1,'  V 
1/  V 

i  i  n  m  »  «  » 

/  '  •  '  9  9 

Vi.ov 

1  1 

9 

9 

VO. 8  V  V  1 

lilt 
9  9 

,*0.9V  V  ' 

RSQ=99.0 
MDP (1 )=1 
ALRP( 1 )=0 . PS 
ALBP (2 ) =D . OS 
IB=7 
IX=10 
IP=20 
N=n 

VALUE=0 . 0 
DIV=1 . 0 
A=0 

NTEST=P 

MTESTOO 


C  DELETE  DATA  PRINTOUT  FILE  IF  FILE  CURRENTLY  IN  EX^ STANCE 

OPEN  1," ASSAY. DATA" ,ATT="OP" ,LEN=] 30 
CLOSE  1 

DELETE  "ASSAY. DATA" 

C  PROGRAM  INTRODUCTION  AND  DESCRIPTION  AS  DISPLAYED  ON  TERMINAL 

TYPE”<NL><NL><NL>  ASSAY  DATA  ANALYSIS  PROGRAM" 

TYPE"<NL>PURPOSE:  FIT  A  LINEAR  OR  POLYNOMIAL  REGRESSION  EOUATION" 
TYPE"  TO  A  STANDARD  CURVE  USED  TO  CALIBRATE  a  CHEMICAL  OR" 

TYPE"  MICROBIOLOGICAL  ASSAY.  A  MINIMUM  OF  FIVE  DATA" 

TYPE"  POINTS  ARE  REQUIRED  FOR  A  LINEAR  REGRESSION  ANALYSIS." 

TYPE"  A  MINIMUM  OF  EIGHT  DATA  POINTS  ARE  RECOMMENDED  FOR" 

TYPE"  A  POLYNOMIAL  REGREWWION  ANALYSIS  (MAXIMUM  DEGREE  OF" 

TYPE"  THE  CALCULATED  POLYNOMIAL  IS  4) .  AFTER  REGRESSION" 

TYPE"  CURVE  HAS  BEEN  CALCULATED,  A  MAXIMUM  OF  SO  UNKNOWNS" 

TYPE"  (ABSORBANCES)  MAY  BF,  ENTERED.  CONCENTRATION  OF  UN-" 

TYPE"  KNOWN  CORRESPONDING  TO  ENTERED  ABSORBANCE  WILL  THEN" 

TYPE"  RE  CALCULATED  USING  REGRESSION  EQUATION.  FOLLOW" 

TYPE"  INSTRUCTIONS  DISPLAYED  ON  TERMINAL  FOR  DATA  ENTRY." 
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TYPE"<NL><NL> 


PAUSE  "PRESS  <RETURN>  TO  CONTINUE" 

TYPE " <N  L  > <N  L>STANDARD  CURVE  ENTRY  ROUTINE: <NL>" 

TYPE"<NL>ENTER  STANDARD  CURVE  DATA  POINTS  AT  THTS  TIME,  IN" 

TYPE" INCREASING  'STANDARD  VOLUME'  ORDER.  WHEN  LAST  DATA  POINT" 

TYPE "HAS  BEEN  ENTERED,  IT  IS  NECESSARY  TO  ENTER  '999'  TO  TERMINATE" 
TYPE " STANDARD  CURVE  ENTRY  ROUTINE. <NL>" 

C  STANDARD  CURVE  DATA  POINT  ENTRY  ROUTINE 

20  N=N+1 

TYPE" - <NL>" 

24  write (in, inno)  n 

WRITE (10, 1005) 

READ  FREE  (11)  VALUE 
DAT  (N ,  1 )  =VALUE 
WPTTE (10,1010) 

READ  FREE  (11)  VALUE 
DAT  (N ,  2 )  =VALUE 

TYPE"<NL>EXAMTNE  ABOVE  ENTRIES.  IF  YOU  WISH  TO  CORRECT" 

TYPE "EITHER  OR  BOTH  VALUE (S),  ENTER  BELOW. <NL>" 

25  TYPE” . ENTER  '1'  TO  CORRECT  DATA  POINT  ENTRY” 

TYPE" . ENTER  ’2’  TO  CONTINUE  WITH  DATA  POINT  ENTRY" 

TYPE" . ENTER  '999*  TO  TERMINATE  DATA  POINT  ENTRY" 

ACCEPT" <NL>ENTER  CHOICE:  ",  I 

IF(I.EQ.l)  GO  TO  24 
IF(I.E0.2)  GO  TO  28 
IF (I. EQ. 999)  GO  TO  30 
GO  TO  25 

28  GO  TO  20 

30  TYPE"<NL><NL><NL><NL>" 

TYPE" ENTER  TYPE  OF  REGRESSION  CURVE  DESIRED  —  LINEAR  OR  POLYNOMIA 
1L<NL>" 

TYPE" . ENTER  '1'  FOR  LINEAR  REGRESSION  CURVE." 

TYPE" . ENTER  '2'  FOR  POLYNOMIAL  REGRESSION  CURVE." 

ACCEPT"<NL>ENTER  CHOICE:  ",  ITEST 

IF(ITEST.EQ.l)  MDP (1 ) -1 
IF(ITEST.EQ. 7)  MDP(1)=4 
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IF(ITEST.NE. l.AND. ITEST.NE.2)  GO  TO  30 


C  DETERMINE  MAXIMUM  VOLUME  (MAXVOL)  AND  MAXIMUM  ABSORBANCE 
C  (MAXABS)  OF  STANDARD  CURVE. 

34  MAXABS=0 . 

MAXVOL=0. 

DO  50  1=1 ,N— 1 

IF(DAT(I+1 ,1) .GT.DAT(I,1))  MAXVOL=DAT(I+l , 1) 

IF(DAT (1+1,2) .GT.DAT (1,2))  MAXABS=DAT(I+1 ,2) 

50  CONTINUE 

C  SET  A=1  IF  'MAXABS'  >1.0 

IF(MAXABS.GT. 1 . )  A=1 

C  DIVIDE  INPUT  STANDARD  VOLUME  BY  'MAXVOL'  TO  SCALE  RANGE  OF 
C  STANDARD  CURVE  OVER  INTERVAL  0.00  TO  1.00.  TRANSFER  STANDARD 
C  CURVE  FROM  ARRAY  'DAT'  TO  ARRAY  'XYW'  FOR  INPUT  TO  SUBROUTINE 
C  'RLFOR* .  SET  *XYW(I,3)'  =  1.0  FOR  UNWEIGHTED  REGRESSION  COEFFICIENTS. 


DO  60  1=1, N 

DAT  (1 , 3)=DAT  (1 , 1)  /MAXVOL 
XYW(I ,1)=DAT(I , 3) 

XYW(I ,2)=DAT(I ,2) 

XYW(I,3)=1.0 
60  CONTINUE 

C  CALCULATE  REGRESSION  EQUATION  USING  IMSL  SUBROUTINE  'RLFOR'. 

70  TYPE " <NL><NL>CALOULATING  POLYNOMIAL  REGRESSION  EQUATION. <ML><NL>" 

CALL  RLFOR  (XYW,IX,N,RSQ,MDP,ALBP,ANOVA,B,IB,PRED,TP,WK,IER) 

C  ABORT  PROGRAM  IF  'IER'  GREATER  THAN  100 

IF(IER.GE.IOO)  STOP  "INADEQUATE  NUMBER  OF  DATA  POINTS" 

C  TRANSFER  REGRESSION  EQUATION  COEFFICIENTS  TO  ARRAY  'COEFF'. 

COEFF ( 1 ) =B ( MDP (2 )+l , 2) 

DO  80  1=1, MDP (2) 

COEFF(I+l)=B(I ,2) 

80  CONTINUE 

C  IF  LINEAR  REGRESSION  EQUATION  CALCULATED  DO  NOT  DETERMINE 
C  MONOTINICITY  OF  EQUATION. 

WRITE (10, 1210)  MDP (2) 
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R5  IFdTEST.En.  l)  GO  TO  100 


C  IF  POLYNOMIAL  REGRESSION  EQUATION  CALCULATED  USE  SUBROUTINE 
C  'MONOT'  TO  DETERMINE  IF  EQUATION  IS  MONOTONIC. 

CALL  MCNOT  (C0EFF,MDP(2) ,XHT ,XLO, IERROR) 

C  IF  REGRESSION  EQUATION  NOT  MONOTONIC  NOTIFY  USER  TO  RECHECK 
C  STANDARD  CURVE  DATA  FOR  CORRECTNESS  OR  SUGGEST  HE  PREPARE  NEW 
C  STANDARD  CURVE. 

IF(IERROR.EQ.O)  GO  TO  100 
NTEST=1 

TYPEM<NL><NL><NL>CAUTION:  CALCULATED  REGRESSION  EQUATION  IS  NOT" 
TYPE"MONOTONIC .  SUGGEST  USER  RE-CHECK  STANDARD  CURVE  DATA  FOR" 
TYPE" ACCURACY.  IF  STANDARD  CURVE  DATA  ENTERED  CORRECTLY  USER" 

TYPE "SHOULD  CONSIDER  PREPARING  NEW  STANDARD  CURVE.  UNKNOWN" 

TYPE "SAMPLE  DATA  WILL  NOT  BE  ACCEPTED  BY  THE  PROGRAM.  STATISTICAL" 
TYPE "DATA  AND  PLOT  OF  REGRESSION  EQUATION  IS  AVAILABLE  TO  USER" 
TYPE" IN  FILE  NAMED  'ASSAY. DATA* .  <NL><NL->" 


PAUSE  "PRESS  <RETURN>  TO  CONTINUE" 

GO  TO  10? 

C  INPUT  AT  TERMINAL  ASSAY  INFORMATION  {NAME,  DATE  OF  ASSAY, TYPE, 

C  INVESTIGATOR  AND  INFORMATION  REQUIRED  FOR  UNKNOWN  CONCENTRATION 
C  CALCULATION . 

100  WRITE (10, 1?S0) 

READ (11,. 1750)  (HD1 (I) ,1=1,4) 

WRITE (10, 1235) 

READ (11, 1250)  (HD2(I) ,T=1,4) 

WRITE  (10, 12.40) 

READdl ,  1250)  (HD3(I)  ,1=1,4) 

WRITE (10,1245) 

READdl, 1250)  (HD4  (I)  ,1=1,4) 

WRITE (10 , 1255) 

WRITE (10,1250) 

READ  FREE  (11)  STAND 
WRITE (10,1255) 

READdl, 1270)  HD 5 
WRITE (10, 1275) 

READ  FREE  (11)  D ILF AC 

IF(DILFAC.GE. 1000.)  GO  TO  130 
GO  TO  135 
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130 


WRITE  (10,12*76) 

WRITE (10,1277) 

READ (11, 1270)  HD6 
DIV=1000. 

135  WRITE (10, 1780) 

READ (11, 1285)  NVOL 

IF (NVOL. LT. 1 .OR.NVOL.GT. 41  GO  TO  135 

DO  140  1=1, NVOL 
WRITE (10, 1290)  I 
READ  FREE  (11)  SVOL 
SVOLUME(I)=SVOL 
140  CONTINUE 

C  INITIALIZE  ELEMENTS  OF  'CONC'  ARRAY  TO  ZERO. 

DO  145  1=1,50 
DO  145  ,1=1,8 
C0NC(I,J)=0. 

145  CONTINUE 

TYPE"<NL><NL>SAMPLE  DATA  ENTRY  ROUTINE<NL>" 


TYPE"  1.)  ON  PROMPT  -  ENTER  SAMPLE  NUMBER." 

TYPE"  2.)  ON  PROMPT  -  ENTER  'ABSORBANCE'  AT  INDICATED  DILUTION" 

TYPE"  VOLUME." 

TYPE"  3.)  IF  YOU  ARE  MISSING  AN  'ABSORBANCE'  AT  INDICATED" 

TYPE”  DILUTION  VOLUME,  ENTER  '0 ' .<NL><NL>" 

KT=0 

150  KT=KT+1 

TYPE" - <NL>" 


155  WRITE (10, 1295)  KT 

WRITE (10, 1300) 

READ  FREE  (11)  SNUM 
SNUMBER (KT ) =SNUW 

DO  180  1=1, NVOL 
WRITE (10, 1305)  I 
READ  FREE  (11)  SABS 
ABSORBANCE (KT , I ) =SARS 
IF (ABSORBANCE (KT, I) .GT.O.)  GO  TO  180 
C0NC(KT,I)=0.0 
GO  TO  180 

160  CONTINUE 

C  USE  SUBROUTINE  'BSECT*  TO  CALCULATE  'STANDARD  ECU I VALENT  VOLUME’ 
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C  FROM  'UNKNOWN'  SAMPLE  ABSORBANCE. 

CALL  BSECT  (SABS, X, COE FF,MDP (2) ,XLO,XHI ,XTOL, IERR) 

C  IF  IERR=f),  SET  CONC  (KT, I)=X.  IF  IERR=1,  SET  C0NC(KT,I)=-111 . 

C  IF  IERR=2 ,  SET  CONC (KT, I )=-222 . 

IF(IERR.EO.O)  CONC (KT, I)=X 
IF(IERR.EQ.l)  CONC(KT,I)=-lll. 

TF(IERR.EO.  2)  CONC (KT, I )=-222. 

180  CONTINUE 

IF(KT. EQ. 50)  GO  TO  184 
GO  TO  188 

184  TYPE"<NL><NL>CAUTION:  YOU  HAVE  ENTERED  ABSORBANCE  DATA  FOR  SO" 

TYPE "SAMPLES.  YOU  MUST  NOW  ENTER  '999'  BELOW  AND  TERMINATE  DATA" 
TYPE "ENTRY.  TO  ENTER  MORE  ABSORBANCE  DATA  USING  SAME  STANDARD" 
TYPE "CURVE  USE  PROGRAM  CONTINUATION  OPTION  TO  BE  PRESENTED" 

TYPE" BELOW.  <NL><NL>" 

IRS  TYPE"<NL>EXAMINE  ABOVE  DATA.  IF  YOU  WISH  TO  CORRECT" 

TYPE "ANY  ENTRY,  ALL  DATA  FOR  SAMPLE  ",SNUM,"  MUST" 

TYPE "BE  RE-ENTERED. <NL>" 

TYPE" . ENTER  '1 '  TO  CORRECT  DATA  ENTRY." 

TYPE" . ENTER  '2'  TO  CONTINUE  WITH  DATA  ENTRY." 

TYPE" . ENTER  '999’  TO  TERMINATE  DATA  ENTRY." 

190  ACCEPT"<NL>ENTER  CHOICE:  ",  I 

IF(I.FQ.l)  GO  TO  155 
IF(I . EO. 2)  GO  TO  150 
IF(I.EQ.999)  GO  TO  205 
GO  TQ  190 

C  CALCULATE  SAMPLE  'UNKNOWN '  CONCENTRATION  FROM  'STANDARD  VOLUME 
C  EQUIVALENT'  OF  'UNKNOWN'  OBTAINED  FROM  'BSECT'  SUBROUTINE. 

205  DO  300  1=1 ,KT 

TEMP1=0 . 

K=0 

DO  240  J=1 ,NVOL 

IF(CONC(T,J) .GT.O.)  GO  TO  210 
GO  TO  240 

210  CONC (I ,J)= (CONC (I , J) *MAXVOL*STAND* ( 1 ,/SVOLUME (J ) ) *DILFAC) 

CONC  (I  ,.7)=CONC(I  ,J)/DTV 
TEMPI =TEMP1+CCNC (I ,J) 

K=K+1 

240  CONTINUE 
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CCNC(I,5)=TEMP1/K 


C  CALCULATE  +/- (MEAN  CONCENTRATION) *0.  IF  ANY  INDIVIDUAL  'UNKNOWN' 

C  CONCENTRATION  EXCEEDS  +/- (MEAN  CONC  +  (MEAN  CONC) *0)  DO  NOT  INCLUDE 
C  THAT  'UNKNOWN'  CONCENTRATION  IN  'AVERAGE  MEAN  ADJUSTED  CONCENTRATION ' . 

MCPLUS=CONC ( I , 5) + (CONC ( I , 5) *0) 

MCMINUS=CONC (I , 5) - (CONC (I , 5) *0) 


TEMP2=0. 

K=0 

DO  270  L=1,NV0L 

IF(CONC(I,L) .GT.O)  GO  TO  260 
GO  TO  270 

260  IF ( (CONC (I  ,L) .GT.MCPLUS)  .OR. 

1  (CONC (I ,  L)  . LT.MCMfNUS)  )  GO  TO  270 

TEMP2=TEMP24C0NC  ( I ,  L) 

K=K+1 

270  CONTINUE 

M=NVOL-K 

IF(K.EO.O)  GO  TO  275 
GO  TO  280 

275  CONC(I,7)=-333. 

GO  TO  290 

280  CONC(I,7)=TEMP2/K 

290  CONC (I ,fi)=M 

C0NC(I,8)=K 

300  CONTINUE 

C  TEST  TO  DETERMINE  IF  'MTEST'  FLAG  HAS  BEEN  CHANGED  TO  '2'.  TF 
C  NOT,  OPEN  FLIE  FOR  INITIAL  DATA  ENTRY.  IF  ' MTEST ’=  2,  BRANCH  TO 
C  LINE  450  TO  ADD  NEW  OUTPUT  DATA  TO  FILE. 

IF (MTEST. EQ. I)  GO  TO  450 

C  OPEN  CHANNEL  M'  FOR  USE  IN  CREATING  DATA  OUTPUT  FILE. 

302  OPEN  1 , " ASSAY . DATA" , ATT= "OP " , LEM=1 30 

C  CREATE  FILE  OF  OUTPUT  DATA  NAMED  ' ASSAY. DATA ' . 

WRITE (1,1015) 

WRITE (1,1018) 

WRITE (1,1020)  MDP(l) 

WRITE (1,1022)  MDP (2) 

WRITE (1,1025) 
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WRITE  (1,1(130) 

WRITE (1,1031) 

WRITE (1,1035) 

WRITE (1,1038) 

WRITE (1,1040)  ANOVA(l) ,ANOVA(4) ,ANOVA(7) ,ANOVA(9) ,ANOVA(10) 
WRITE (1,1043)  ANOVA (2) ,ANOVA(5) ,ANOVA(R) 

WRITE (1,1045)  ANOVA ( 3 ) , ANOVA ( 5 ) 

WRITE (1,1030) 

WRITE (1 , 1047) 

WRITE (1 , 1050) 

WRITE (1,1 053) 

WRITE (1,1055) 

WRITE (1,1057) 

WRITE (1 ,1058) 


no  400  J=1,MDP(2) 

WRITE (1 , 1050)  J,  (B(J,I) ,1=1,8) 
400  CONTINUE 


WRITE (1 , 1053)  (B(MDP(2)+1,I) ,1=1,8) 


WRITE (1,1050) 

WRITE (1,10*5)  ANOVA (11) 
WRITE (1, 1J57)  ANOVA (12) 
WRITE (1,1070)  IER 
WRITE (1,1085) 

WRITE (1,1090) 

WRITE (1,1093) 

WRITE  (1,1.030) 

WRITE (1 ,1095) 

WRITE (I  ,1097) 

WRITE (1,1 098 ) 

WRTTE (1 ,11  00) 

WRITE (1,11 03) 


DO  420  1=1, N 

RES=DAT (I , 2)-PRED (1,1) 

WRITE  (1. ,  11  05)  I  ,DAT(T  ,  1)  ,DAT(I,3)  ,DAT(I,2)  ,PRED(I,1)  ,RFS 
420  CONTINUE 


WRITE '1, 1030) 

WRITE (1,1107) 

WRITE (1,1) 10) 

C  CALL  PLOT1  AND  PLOT?  SUBROUTINES  TO  CREATE  LOW  RESOLUTION  GRAPH 
C  OF  INPUT  STD  CURVE  DATA  POINTS  AND  CALCULATED  REGRESSION  EQUATION. 

CALL  PLOT1  (GRAPH) 

CALL  PLOT 2  (GRAPH,N,DAT,MDP(2) ,COEFF,A) 
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WRITE  (1 , 1085) 
WRITE (1,1430) 
WRITE (1,1432) 
WRITE (1,1434) 


.1=50 

1=1 


430  WRITE (1 , 1435)  LABEL1 (I ) ,TABEL?  (I ) , (GRAPH (J,K) ,K=0 , 100) 

J=vJ-l 

1=1+1 

IF(.l.GE.O)  GO  TO  430 
CONTINUE 

IF(A.EQ.O)  GO  TO  435 
GO  TO  440 

435  WRITE (1 , 1440) 

GO  TO  445 

440  WRITE (1,1442) 

445  WRITE (1, 1445) 

C  IF  REGRESSION  EQUATION  NOT  NONOTONIC  DO  NOT  ADD  UNKNOWN  CONCENTRATION 
C  DATA  TO  OUTPUT  DATA  FILE. 

IF (NTEST. EQ. 1)  GO  TO  300 

450  WRITE (1,1310) 

WRITE (1,1315) 

WRITE (1 , 1320)  (HD1(I) ,1=1,4) , STAND, HD5 
WRITE (1 , 1325)  (HD2(I) ,1=1,4) ,DILFAC 
WRITE (1 , 1330)  (HD3(I),I=1,4) ,NVOL 

IF(DILFAC.GT. 1000 . )  GO  TO  304 
GO  TO  300 

304  WRITE (1 , 1335)  (HD4 (I) ,1=1 ,4) ,HP5 

GO  TO  308 

305  WRITE (1,1335)  (HD4 (I ) , 1=1 , 4) ,HD5 

308  DO  110  1=1 ,NVOL 

WRITE (1,1 340)  I,SV0LUME(T) 

310  CONTINUE 

WRITE (1,1 34 5) 

WRITE (1,1 350) 

WRITE (1,1 355) 

WRITE (1,1 350) 

WRITE (1,1 355) 

WRITE (1,1 370)  0 
WRITE (1,1 375) 
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no  140  T  =  |  ,trr 

'•••PTTK(1  ,  n«0)  PMUMBFR  (T )  , (ABSORBANCE (I ,0)  ,,7=1  ,4)  , 

1  (CONC(I,K) ,K=1,R) 

340  r<  INI  IE 

WRITE  (!  ,  I  350) 

WRITE  (]  ,1W) 

WRITF  (1,13^0) 

WRITE (1,1395) 

WRITE (1,1400) 

WRITE (1,1415) 

WRITE (1 ,1420)  n 

350  TYPF”<NL><NL><NL><N[>" 

TYPE "DO  YOW  WISH  TO  ENTER  MORE  DATA  USING  SAME  STANDARD  CURVE?< 
1NL>" 

TYPE" . ENTER  '  1 '  TO  CONTINUE  DATA  ENTRY" 

TYPE" . ENTER  '2'  TO  EXIT  PROGRAM" 

ACCEPT" <NL>ENTER  CHOICE:  " ,  MTEST 
TYPE"<NL><NL><NL><NL>" 

IF (MTEST. E0.1)  GO  TO  100 
IF (MTEST. E0- 2)  GO  TO  ^60 
IF(MTEST.NE.1.AND.MTEST.NF,.7)  GO  TO  350 

360  CLOSE  1 

C  TERMINAL  DATA  TNPIJT  AND  STATISTICAL  OUTPUT  FORMAT  STATEMENTS. 

1000  FORMAT ( / , "DATA  POINT  MO.:" ,13,/) 

1005  FORMAT (Z, "ENTER  VOLUME  OF  STANDARD  ('IND'  VARIABLE):  ") 

1010  FORMAT (Z, "ENTER  MEASURED  ABSORBANCE  FOR  ABO’Tl  VOLUME (*DEP’  VARIAPL 
C1E) :  ") 

1015  FORMAT (IX,"***  STATISTICAL  DATA  — >  CALCULATED  POLYNOMIAL  REGPFSSI 
ION  EQUATION  ***") 

101ft  FORMAT (IX," - 

1 - ",//) 

1020  FORMAT (IX, "REQUESTED  DEGREE  OF  REGRESSION  EQUATION:  ”,T1) 

1022  FORMAT (IX, "DEGREE  OF  CALCULATED  REGRESSION  EQUATION:  "  ,11,/) 

1025  FORMAT (IX,"**  ANALYSIS  OF  VARIANCE  **") 

1030  FORMAT (IX," - 

1 - ”) 

1033  F0RMAT(1X,”*",7X, "SOURCE  OF" ,PX, "DEGREES  OF", 3X, "SUMS  OF", 5X, "MEAN 

1 ",SX, "COMPUTED" ,6X,"TAIL  AREA" ,6X,"*") 

1035  FORMAT (IX , "* " , 2X , "VAF IATION" ,9X , " FREEDOM" , 5X , " SQUARES" , 4X , "SQUARES 
1" , 4X , " F-VALUE" , 5X , " F-DISTR IBUT ION" , 3X , "* ”) 

1038  FORMAT  (IX,  3X," - ",8X," - ",3X," - ",4X," - " 

l,4X," - ",6X," - ",/) 

1040  FORMAT  (IX,  3X,  "REGRESSION", 9X,F4.0,5X,F9. 5, 7X,F9. 5, 3X,FR.  3,4X,F14.9 


1) 

1041  FORMAT (IX,  3X,  "RESIDUAL"  ,  11X, F4. 0, 5X, F9. 5, 2X,F9.  5) 

1045  FORMAT ( IX, 3X, "CORRECTED  TOTAL" ,4X,F4.0,5X,F9.5) 

1047  FORMAT (//////, IX,"**  REGRESSION  COEFFICIENTS  **“) 

1050  FORMAT (IX," - 

1053  FORMAT  (IX ,  "* " , 50X , " LOWER"  ,  7X ,  "UPPER"  , 5X ,  "ADJUSTED"  , 3X ,  "PARTIAL"  ,  17 
IX , "*") 

1055  FORMAT (IX , , "DEGREE" , 5X , "VARIABLE"  , 4X , "REGRESSION" , 4X , " STANDARD" 
1 , 3X, "CONFIDENCE" , 2X , "CONFIDENCE" , 3X , "SUMS  OF" , 4X , "F-TEST" , 5X , "TAIL 
1  AREA" ,4X, "*") 

1057  F0RMAT(1X,"*" ,"OF  TERM" ,5X, "MEANS" ,4X, "COEFFICIENTS" ,4X, "ERROR" ,7X 
1 , "LIMIT" , 7X , "LIMIT" ,6X , "SQUARES" , 4X , "VALUES" , 3X , " F-DTSTRTBUTTON" , 1 
IX,"*") 

1058  FORMAT  (2.X , " - " , 4X , " - "  , 3X , " - " ,  3X , " - " , 

13X," - ",2X," - "  ,  3X , " - "  ,  3X , " - ",2X," - 

1 - ",/) 

1050  FORMAT(lX,lX,"X**",Il,7X,F8.5,2X,F13.8,3X,FR.5,3X,F8.4,4X,F8.4,4X, 

1F9.5,2X,F8.4,2X,F14.9) 

1083  FORMAT (IX,  IX, 11  INTERCEPT: 11 ,1X,F8.5,7X,F13.8,3X,F8.5,3X,F8.^,4X,F8.4 
1,4X,F9.5,2X,F8.4,2X,F14.9) 

1085  FORMAT (//, IX ."PERCENTAGE  OF  VARIATION  EXPLAINED  BY  ESTIMATED  MODEL 
1:  " ,F10. 6) 

1087  FORMAT (IX, "STANDARD  DEVIATION  OF  RESIDUALS:  " ,F10.S,/) 
lO^O  FORMAT ( 1 X,"PO LYNOM I AL  CALCULATION  SUBROUTINE  ERROP  PARAMETER:  ",T1 
1) 

1085  FORMAT ( 1 H 1 ) 

1090  FORMAT (IX,"***  INPUT  DATA  AND  PREDICTED  VALUES  FROM  MODEL  ***") 

1093  FORMAT  (IX," - ",/// 

1) 

1 095  FORMAT ( IX , " * " , 17X , "-X- " ,9X , "-X-" , 51X , " * " ) 

1 097  FORMAT (IX , "* " , 16X , " INPUT" , 7X , "SCALED" ,9X , "-Y-" , 1 IX , "-Y-" , 23X , "* ") 

1098  FORMAT (IX,"*" ,2X,"DATA  POINT" ,3X," STANDARD” ,4X, "STANDARD" ,7X,"INPU 
IT" ,7X, "PREDICTED" ,20X,"*") 

1100  FORMAT  ( 1  X ,  "* "  , 4X ,  "NUMBER"  ,6X  ,  ''VOLUME"  ,4X,"EQ.  VOLUME”  , 4X , " ABSORRAN 
ICE" ,4X, "ABSORBANCE" ,4X , "RESIDUAL  VALUE" , IX,"*") 

1103  FORMAT  (IX,  3X, " - "  , 3X , " - "  , 3X , " - " , 4X , " - 

1 _ "  ,4X,  " _ "  ,4X," _ "  ,/) 

1105  F0RMAT(1X,7X,I7,7X,F7.4,SX,F‘7.4,RX,F5.4,8X,F8.4,9X,F7.5) 

1107  FORMAT (/, IX, 7X,”X  — >  INDEPENDENT  VARIABLE") 

1110  FORMAT (IX, 7X,"Y  — >  DEPENDENT  VARIABLE") 

1.210  FORMAT  ("POLYNOMIAL  REGRESSION  EQUATION  OF  DEGREE  ",Il,"  HAS  BEEN  C 
1ALCULATED." ,//) 

C  FORMAT  STATEMENTS  FOR  UNKNOWN  SAMPLE  DATA  OUTPUT  AND  INPUT. 

1230  FORMAT (Z," ENTER  TYPE  OF  ASSAY:  ") 

1235  FORMAT (Z," ENTER  CATE  OF  ASSAY:  ”) 

1240  FORMAT (Z, "ENTER  TYPE  OF  SAMPLE:  ") 

1245  FORMAT (Z, "ENTER  NAME  OF  INVESTIGATOR:  ") 
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1  250 
1255 
1260 

1265 

1270 

1275 

1276 


1277 

1280 


l?85 
1790 
1295 
1300 
1  305 
1310 
1315 

1320 


1325 


1330 

1335 

i  i/in 
1345 
1  3S(J 


1  355 

i  JGp 


1765 

1370 


1375 


1380 


FORMAT (4 A8) 

FORMAT ( / , " STANDARD  CONCENTRATION  DATA:  ") 

FORMAT {" . ENTER  NUMERICAL  PORTION  OF  STANDARD  CONCENTRATION:  ", 

12) 

FORMAT (" . ENTER  UNITS  OF  STANDARD  (EXAMPLE:  GM/ML) :  "  ,Z) 

FORMAT (AS) 

FORMAT (//'ENTER  INITIAL  DILUTION  FACTOR  OF  SAMPLE  MATERIAL:  ",Z) 
FORMAT (//'CAUTION:  DILUTION  FACTOR  GREATER  THAN  1000" ,/,9X, "YOU  MU 
1ST  INCREASE  UNITS  OF  STANDARD  BY  A  FACTOR  OF  1000 . " ,//,9X,"IF  UNIT 
IS  OF  STANDARD  WERE  'NG/ML' ,  CHANGE  TO  'UG/ML'  ’’,/,9X,"IF  UNITS  WER 
IE  'UG/ML',  CHANGE  TO  'MG/ML'  ",/,9X,"IF  UNITS  WFRE  'MG/ML',  CHANGE 
1  TO  'GM/ML'  ",/) 

FORMAT (Z," . ENTER  NEW  UNITS  FOR  STANDARD:  ’’) 

FORMAT  (//'ENTER  NUMBER  OF  SAMPLE  DILUTIONS  PER  SAMPLE  (MINIMUM  1, 

1  MAXIMUM  4):  ",Z) 

FORMAT (II) 

FORMAT ( " . ENTER  'MLS  OF  SAMPLE'  AT  DILUTION  ",I1,"  :  ",Z) 

FORMAT  (//'DATA  ENTRY:  ",I3,/) 

FORMAT ("ENTER  SAMPLE  NUMBER:  H,Z) 

FORMAT ( " . ENTER  ABSORBANCE  AT  SAMPLE  DILUTION  ",I1,":  ",Z) 

FORMAT (1H1 ,"***ASSAY  REGRESSION  PROGRAM  -  SAMPLE  nATA  OUTPUT***") 

FORMAT  (IX," - ",/ 

1/) 

FORMAT (IX /'TYPE  OF  ASSAY:" ,8X,4A8,9X, "ASSAY  PARAMETERS  ->  STANDARD 
1  CONCENTRATION :",1X,F7. ?, IX, A6) 

FORMAT (IX, "DATE  OF  ASSAY:" ,8X,4A8,9X," -  SAMPLE  D 

1 1 LUTION  FACTOR ; " , IX ,  F7 . 2 ) 

FORMAT (IX, "SAMPLE  TYPE:" ,10X,4A8,29X,"NO.  OF  SAMPLF  DILUTIONS: " ,1X 
1,11) 

FORMAT (IX, "NAME  OF  INVESTIGATOR : " , IX , 4A8 , 29X , " SAMPLE  CONCENTRATION 
1  UNITS:  " ,A6,/) 

FORMAT (IX, R3X," SAMPLE  DILUTION  ",I1,"  =  ",F5.3,"  ML  ") 

FORMAT (/) 

FORMAT (IX," - 


1-") 

FORMAT  (IX , "* "  ,  1 .1 X ,  "* " , 31X , "* " , 31 X , "*"  ,8X , "*" , 5X ,  "ADTUSTFD  CONCENTR 
1 AT TOM" ,8X , "*") 

FORMAT (IX,"*" ,1 IX, "*" ,6X,"A  BSORBANC  E" ,SX,"*" ,3X,"C  0  N  C 
1  ENTRATIO  N" , 3X,"*" ,3X,"MEAN" ,2X , "*" ,2X, "NO.OF  DTL  VOLS",7X 
1 , "* " , IX , "ADJUSTED" , IX , "* " , IX , "NO. " , IX , "* ") 

FORMAT (IX,"*", IX, "SAMPLE  NO.*" ,31X,"*"  ,31X,"*" ,2X,"CONC" ,2X,"*",4X 
1, "EXCEEDING", 5X,"*",3X, "MEAN", 3X/'*",  IX, "IN", 2X,"*") 

F0RMAT(1X,"*",11X,"*",1X,"DIL  1  DIL  2  DIL  3  DIL  4", IX,"*", IX 
I, "DILI  DIL  7  DIL  3  DIL  4",1X,"*",2X,"  (MC)  ",2X,"*",.1X,"+/-  (" 
.1 ,  F3 . 2,  "*MC  +  MC)  ", IX, "*",3X,"CONC",3X,"*", IX, "AVG. ","*") 

FORMAT  (IX,  IX," -  -  -  -  -  - 


1  — ",/) 

FORMAT (.1 X,  3X, I4,7X , F5. 3, 3X,  F5. 3, 3X ,F5.  3, 3X,F5. 3, 2X,F6. 1  , 2X,F6. 1 , ?X 
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1 ,  F6 . 1 , 2X ,  F5 . 1 , 3X ,  F5 . 1 ,  IX , "  f  “  ,  RX ,  F? .  0 , 11 X ,  F5 . 1 , 4X ,  F2 . 0 , 2X , "  1  ") 

1385  FORMAT (//) 

1390  FORMAT (IX, "ERROR  CODES  -  CONCENTRATION'  CALCULATION  ROUTINE") 

1395  FORMAT  (IX," - ",/) 

1400  FORMAT (IX,"  1.)  -111.  IN  CONCENTRATION  COLUMNS  INDICATES  THAT 

1  SAMPLE  ABSORBANCES  NOT  IN  RANGE  OF  STANDAPD  CURVE.",/) 

1415  FORMAT (IX,"  2.)  -333.  IN  ADJUSTED  CONCENTRATION  COLUMNS  INDTCA 

1TES  TOAT  THE  CONCENTRATION  OF  UNKNOWN  AT  ALL  DILUTIONS") 

1420  FORMAT (IX , "  EXCEEDED  ->  +/- (" ,F3. 2,"*MC  +  MC) .") 

1430  FORMAT  (IX, 18X,"***  CALCULATED  REGRESSION  EQUATION  -  LOW  RESOLUTION 
1  PLOT  -  [X-AXIS  VERTICAL;  Y-AXIS  HORIZONTAL!  ***”,/) 

1432  FORMAT ( IX, 82X," SYMBOL  FOR  INPUT  DATA  POINTS:  ->  0") 

1434  FORMAT ( IX, 82X," SYMBOL  FOR  REGRESSION  EQUATION:  ->  «",/) 

1435  FORMAT (IX, 1 OX, A1 ,2X,A3,1X,101A1) 


1440 

FORMAT  (IX,  16X,  "0.0 

0.1 

0.2 

0.3  0.4 

n 

1.5  0.5 

0.7 

O.R 

0.9 

1.0",/) 

1442 

FORMAT ( IX, 15X, "0.0 

0.2 

0.4 

0.5  O.R 

1 

1.0  1.2 

1  .4 

1.5 

1.8 

2.0",/) 

1445 

FORMAT (IX, 57X, "A  B 

S  0 

R  B  A 

N  C 

E"  ,/) 

9000  STOP  "OUTPUT  DATA  FILE  'ASSAY. DATA'  CREATED.  END  OF  PROGRAM" 
9999  END 
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,J  *#######*##*#**##*##*####**#*#***#*#********##*«•##***#*********♦♦**#*# 


r  *  * 

C  *  SUBROUTINE  ’  BSECT '  -  CALLINC  PROGRAM  NAME:  ’ASSAY’  * 

C  *  * 

C  *  WRITTEN  MAY  1981  -  J.J.KNUDSEN,  LAIR,  DIV.  OF  RESEARCH  SUPPORT  * 

C  *  ANALYTICAL  CHEMISTRY  GROUP.  * 

j  *  * 

C  *  LANCUACE  -  DATA  GENERAL  FORTRAN  V  * 

1  * 

'  *  PURPOSE  OF  SUBROUTINE  * 

r>  •* _ * 


C  *  POLY (X )  IS  A  POLYNOMIAL  IN  X  OF  DEGREE  ’NDEG  ’ .  * 

C  *  * 


C  *  POLY ( X )  =  COEFF( ’ZERO’)  +  SUM  COEFF(l )*X**I  * 
C  *  * 
C  *  WHERE:  I  M  TO  NDEG  * 

C  *  * 

C  *  POLY (X )  IS  DEFINED  FOR  ’X’  BETWEEN  ’XLO’  AND  ’XHI*.  GIVEN  A  * 
C  *  PARTICULAR  NUMBER  ’Y’,  THIS  SUBROUTINE  FINDS  THE  VALUE  OF  ’X’  * 
C  *  FOR  WHICH  ’ Y ’ =POLY (X ) .  IN  OTHER  WORDS  IT  NUMERICALLY  SOLVES  THE  * 
C  *  EQUATION  * 
C  *  * 
C  *  F ( X )  =  POLY (X )  -  ’Y’  =  0  * 
C  *  * 
C  *  FOR  A  UNIQUE  SOLUTION  TO  EXIST,  POLY(X)  MUST  BE  MONOTONIC  AND  ’Y’  * 
C  *  MUST  LIE  BETWEEN  POLY (’XLO’)  AND  POLY(’XHI’).  * 


C  *  DESCRIPTION  OF  ALGORITHM  * 

0  * . - .  * 

C  *  THE  ROOT  OF  THE  POLYNOMIAL  LIES  ON  THE  INTERVAL  ( A, B) * (XLO ,XHI ) .  * 

C  *  THIS  MEANS  THAT  CHANGES  SIGN  ONCE  ON  THE  INTERVAL,  SO  THAT  * 
S  *  F(A)*F(B)  <  0.  PICK  A  POINT  MIDWAY  BETWEEN  A  AND  B  (CALL  IT  X).  * 

C  *  IF  F( A ) *F(X )  <  0,  THE  ROOT  IS  BETWEEN  A  AND  X;  B  IS  THEN  ASSIGNED  * 
0  *  THE  VALUE  OF  X.  IF  F(A)*F(X)  >  0,  THE  ROOT  IS  BETWEEN  X  AND  B;  * 
C  *  A  IS  THEN  ASSIGNED  THE  VALUE  OP  X.  AT  THIS  POINT  WE  KNOW  THAT  THE  * 
S  *  ROOT  IS  BETWEEN  A  AND  B,  BUT  THE  INTERVAL  OF  UNCERTAINTY,  (B-A),  * 

C  *  IS  HALF  AS  LARGE  AS  IT  WAS  INITIALLY.  THE  BEST  GUESS  FOR  THE  ROOT  * 
C  *  IS  (A+B)/2,  AND  THE  MAXIMUM  ERROR  IN  THIS  ESTIMATE  IS  IS  (B-A)/2.  * 

C  *  KEEP  REPEATING  THE  ABOVE  UNTIL  THE  ERROR  IS  SMALLER  THAN  ’XTOL*,  * 
C  *  WHICH  HAS  BEEN  SUPPLIED  TO  THE  SUBROUTINE  AS  A  FORMAL  PARAMETER.  * 
C  *  * 

C  *  IN  SUMMARY,  EACH  ITERATION  HALFS  THE  MAXIMUM  ERROR  IN  THE  ESTIMATE  * 
C  *  OF  THE  ROOT.  THIS  ALGORITHM  DOES  NOT  CONVERGE  AS  FAST  AS  SOME  * 
C  *  MORE  SOPHISTICATED  METHODS  (E.G.,  NEWTON’S  METHOD  OR  MODIFIED  * 

C  *  REGIJLA  FALSI),  BUT  IT  IS  GUARANTEED  TO  FIND  THE  ROOT.  OTHER  ALGOR-  * 
C  *  1THMS  ARE  SOMETIMES  TEMPERAMENTAL.  * 

C  *  * 

C  *  SUBROUTINE  FORMAL  PARAMETERS  * 

C  * .  * 
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C  *  INTEGER  * 

C  * .  * 

C  *  IERK  -  SUBROUTINE  ERROR  PARAMETER  * 

C  *  1.)  SET  =  0  IF  NO  ERRORS  ENCOUNTERED  * 

C  *  2.)  SET  =  1  IF  NO  SOLUTION  EXISTS  TO  EQUATION  OVER  * 

C  *  INTERVAL  'XLO'*  0.0  TO  'XHI*-  1.0  * 

C  *  3. )  SET  =  2  IF  NO  SOLUTION  FOUND  TO  EQUATION  IN  MAXIMUM  * 

C  *  NUMBER  OF  ITERATIONS  ( 'MAXITR *  )  * 

C  *  NDEG  -  DEGREE  OF  REGRESSION  EQUATION  GENERATED  BY  CALLING  * 

C  *  PROGRAM  'ASSAY'  * 

C  *  * 

C  *  DOUBLE  PRECISION  * 

C  * . * 

C  *  Y  -  INPUT  'UNKNOWN'  SAMPLE  ABSORBANCE  * 

C  *  X  -  OUTPUT  'STANDARD  VOLUME  EQUIVALENT'  (ESTIMATE  OF  ROOT  OF  * 

C  *  EQUATION  AT  INPUT  'UNKNOWN'  ABSORBANCE).  * 

*  COEFF  -  VECTOR  OF  LENGTH  6  CONTAINING  REGRESSION  COEFFICIENTS  * 

*  XLO  -  LOWER  LIMIT  OF  'STANDARD  VOLUME  EQUIVALENT'  * 

C  *  SET  =  TO  0.0  * 

C  *  Xlil  -  UPPER  LIMIT  OF  'STANDARD  VOLUME  EQUIVALENT'  * 

C  *  SET  =  TO  1.0  * 

C  *  XTOL  -  TOLERANCE  VALUE  FOR  MAXIMUM  ALLOWED  ERROR  IN  CALCULATED  * 

C  *  ROOT.  SET  *  TO  1.0E-P  * 

C  *  * 

C  *  SUBROUTINE  SYMBOL  DEFINITIONS  (EXCLUDING  FORMAL  PARAMETERS)  * 

C  * . * 

C  *  INTEGER  VARIABLES  * 

C  * _ * 

C  *  NDEGP1  -  DEGREE  OF  REGRESSION  EQUATION  +  1  (LOOP  CONTROL)  * 

C  *  MAXITR  -  MAXIMUM  NUMBER  OF  ITERATIONS  ALLOWED  TO  SOLVE  EQUATION  * 

C  #  WITHIN  GIVEN  ERROR  TOLERANCE  * 

C  *  ITR  -  VARIABLE  CONTAINING  COUNT  OF  NUMBER  OF  ITERATIONS  COMPLETED  * 
C  *  BY  ROUTINE  * 

C  *  * 

C  *  DOUBLE  PRECISION  VARIABLES  * 

0* _ _ _  * 

2  *  XERROR  -  VARIABLE  SET  =  TO  XERROK/2  DURING  EACH  ITERATION  OF  * 

C  *  SOLUTION  ROUTINE.  INITIAL  VALUE  =  ARS( ' XHI ' - ' XLO ' .  * 

C  *  A  -  TEMP  STORAGE  VARIABLE.  INITIAL  VALUE  -  'XHI'  * 

C  *  B  -  TEMP  STORAGE  VARIABLE.  INITIAL  VALUE  *  ' XLO '  * 

C  *  POLYA  -  VARIABLE  SET  =  TO  SOLUTION  OF  POLYNOMIAL  AT  'XLO'  * 

C  *  POLYB  -  VARIABLE  SET  =  TO  SOLUTION  OF  POLYNOMIAL  AT  'XHI'  * 

C  *  POLYX  -  VARIABLE  SET  -  TO  SOLUTION  OF  POLYNOMIAL  AT  ROOT  EST  'X'  * 

C  *  FA  -  VARIABLE  SET  =  TO  'POLYA  -  Y’  * 

C  *  FB  -  VARIABLE  SET  *  TO  'POLYP  -  Y’  * 

C  *  FX  -  VARIABLE  SET  *  TO  'POLYX  -  Y’  * 

C  *  * 

n  ********************************************************************** 

SUBROUTINE  B3ECT  (Y*X# COEFF, NDEG  #XLO , XHI t XTOL, I ERR ) 


C  SUBROUTINE  DATA  TYPE  STATEMENTS 

INTEGER  NDEG.NDEGP1 , IERR .MAXITR , ITR 

DOUBLE  PRECISION  Y , X , COEFF ( 6 ) , XLO , XHI , XTOL , XERROR , A , B , FA , FB , 
1  FX , POLY  A , POLYB , POLYX 

C  SUBROUTINE  SYMBOL  INITIALIZATION  STATEMENTS 

ITR-0 

IERR-0 

MAXITR=20 

A-XLO 

B-XHI 

XERROK*=DABS(B-A) 

NDEGP1 'NDEG+1 

C  EVALUATE  POLYNOMIAL  GENERATED  BY  'ASSAY'  USING  HORNERS  METHOD 
C  AT  'XLO'  AND  'XHI'. 

POLYA *0. 

POLYB-O. 

DO  20  I-1.NDEG 

P0LYA«(P0LYA+C0EFF(NDEGP1-I+1 ))*A 
POLYB*=  (P0LYB+C0EFF(NDEGP1  -1  +  1  )  )*B 
20  CONTINUE 

POLY  A  =  POLY  A  +COEFF ( 1 ) 

P0LYB-P0LYB+C0EFF(1 ) 

FA-POLYA-Y 

FB=POLYB-Y 

C  CHECK  FOR  THE  EXISTANCE  OF  A  SOLUTION  BETWEEN  'XHI'  AND 
C  'XLO'.  IF  NO  SOLUTION  EXISTS  SET  'IERR'  TO  '1*  AND  RETURN 
C  TO  'ASSAY'. 

IF(FA*FB. LE.O. )  GO  TO  30 
IERR=1 
RETURN 

30  CONTINUE 

C  CALCULATE  NEW  ESTIMATE  OF  ROOT 

X»(A+B)/2. 

ITR-ITR+1 

C  CHECK  TO  SEE  IF  MAXIMUM  NUMBER  OF  ITERATIONS  EXCEEDED. 
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IF(ITR.LE.MAXITR)  GO  TO  40 
IERR«2 
RETURN 

40  CONTINUE 

C  DIVIDE  MAXIMUM  ERROR  BY  'TWO*. 

XERR0R=XERR0R/2. 

C  CHECK  TO  SEE  IF  MAXIMUM  ERROR  IS  WITHIN  TOLERANCE. 

I F ( XERROR . LT . XTOL )  RETURN 

C  EVALUATE  POLYNOMIAL  AT  NEW  VALUE  OF  ’X'. 

P0LYX=0. 

DO  50  1=1 ,  NDEG 

P0LYX=(P0LYX-C0EFF(NDEGP1 -1+1 ) )»X 
50  CONTINUE 

POLYX=POLYX+CQEFF( 1 ) 

FX=POLYX-Y 

IF(FA*FX. LE.O. )  GO  TO  60 
A=X 
FA=FX 
GO  TO  30 

60  CONTINUE 


B=X 

GO  TO  30 


END 
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Q  ***#*#***##****###**#*******#*******#****#******#***##**#*#*********** 

c  *  * 

C  *  SUBROUTINE  'MONOT'  -  CALLING  PROGRAM  NAME:  'ASSAY'  * 

C  *  * 

C  *  WRITTEN  MAY  1981  -  J.J.KNUDSEN,  LAIR,  DIV.OF  RESEARCH  SUPPORT,  * 

C  *  ANALYTICAL  CHEMISTRY  GROUP.  * 

C  *  * 

C  *  LANGUAGE  -  DATA  GENERAL  FORTRAN  V  * 

C  *  * 

C  *  PURPOSE  OF  SUBROUTINE  * 

C  * .  * 

C  *  THIS  SUBROUTINE  CHECKS  THAT  THE  REGRESSION  POLYNOMIAL  IS  MONOTONIC  * 

C  #  ON  (XLO.XHI)  BY  MAKING  SURE  THAT  THE  DERIVATIVE  DOES  NOT  CHANGE  * 

C  *  SIGN  FOR  100  POINTS  SPACED  EVENLY  ON  THE  INTERVAL.  THIS  SUB-  * 

C  *  ROUTINE  NOT  USED  IF  GENERATED  REGRESSION  EQUATION  IS  LINEAR.  * 

C  *  * 

C  *  SUBROUTINE  FORMAL  PARAMETERS  * 

C  * -  * 

C  *  INTEGER  * 

C  * -  * 

C  *  NDEG  -  DEGREE  OF  REGRESSION  EQUATION  GENERATED  BY  CALLING  * 

C  *  PROGRAM  * 

C  *  IERROR  -  SUBROUTINE  ERROR  PARAMETER  * 

C  *  SET  =■  TO  0  IF  1ST  DERIVATIVE  OF  REGRESSION  EQUATION  * 

C  *  IS  MONOTONIC  * 

C  *  SET  *  TO  1  IF  1ST  DERIVATIVE  OF  REGRESSION  EQUATION  * 

C  *  IS  NOT  MONOTONIC  OVER  RANGE  'X’-O.O  TO  1.0  * 

C  *  * 

C  *  DOUBLE  PRECISION  * 

C  * -  * 

C  *  COEFF  -  VECTOR  OF  LENGTH  6  CONTAINING  REGRESSION  COEFFICIENTS  * 

*  XLO  -  LOWER  LIMIT  OF  'STANDARD  VOLUME  EQUIVALENT’  * 

*  XHI  -  UPPER  LIMIT  OF  'STANDARD  VOLUME  EQUIVALENT’  * 

C  *  * 

C  *  SUBROUTINE  SYMBOL  DEFINITIONS  (EXCLUDING  FORMAL  PARAMETERS)  * 

C  * -  * 

C  *  INTEGER  VARIABLES  * 

C  * -  * 

C  *  NDEGM1  -  DEGREE  OF  REGRESSION  EQUATION  -  1  (LOOP  CONTROL)  * 

C  *  * 

C  *  DOUBLE  PRECISION  VARIABLES  r  * 

C  * .  * 

C  *  DCOEFF  -  VECTOR  OF  LENGTH  6  CONTAINING  COEFFICIENTS  OF  1ST  * 

C  *  DERIVATIVE  OF  REGRESSION  EQUATION  * 

C  *  X  -  VARIABLE  SET  -  TO  VALUE  BETWEEN  0.0  AND  1.0  AT  100  EQUALLY  * 

C  *  SPACED  POINTS.  1ST  DERIVATIVE  EVALUATED  AT  EACH  POINT  * 

C  *  DRVX  -  VALUE  OF  DERIVATIVE  AT  'X'  * 

C  *  DRVXLO  -  VALUE  OF  DERIVATIVE  AT  'XLO'  * 

C  *  TEST  -  VARIABLE  SET  -  TO  'DRVX*DRVXLO ’  TO  CHECK  IF  SIGN  OF  * 
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c  * 
c  * 

C  *■ 


DERIVATIVE  HAS  CHANGED 


* 

****************** 


SUBROUTINE  MONOT  ( COEFF , NDEG ,  XHI , XLO , IERROR ) 

C  SUBROUTINE  DATA  TYPE  STATEMENTS 
INTEGER  NDEG, NDEGM1 .IERROR 

DOUBLE  PRECISION  COEFF(6),DCOEFF(6),XLO,XHI,X,DRVX, 

1  DRVXLO .TEST 

C  CALCULATE  COEFFICIENTS  OF  1ST  DERIVATIVE  OF  POLYNOMIAL 
C  REGRESSION  EQUATION  GENERATED  BY  'ASSAY' 

DO  10  1=1 , NDEG 

DCOEFF( I )=COEFF( 1+1 )*DBLE ( FLOAT ( I ) ) 

10  CONTINUE 

C  CALCULATE  VALUE  OF  1ST  DERIVATIVE  OF  REGRESSION  EQUATION  AT 
C  'XLO'  USING  HORNER'S  METHOD. 

NDEGM1 =NDEG-1 
DRVXL0=0. 

DO  20  1=1 , NDEGM1 

DRVXL0*(DRVXL0*DC0EFF(NDEG-I+1 ) )*XLO 
20  CONTINUE 

DRVXLO=DRVXLO  +  DCOEFF( 1  ) 

C  CALCULATE  VALUE  OF  1ST  DERIVATIVE  OVER  RANGE  'XLO’  TO  'XHI'  AT 
C  100  EQUALLY  SPACED  POINTS.  AFTER  DERIVATIVE  FOR  EACH  POINT 
C  CALCULATED,  TEST  TO  DETERMINE  IF  DERIVATIVE  IS  POSITIVE 
C  IF  POSITIVE  CALCULATE  DERIVATIVE  AT  NEXT  POINT. 

C  IF  NEGATIVE  SET  'IERROR'  TO  1  AND  RETURN  TO  CALLING  PROGRAM. 

DO  40  1=1,100 

X=XLO *0.01 *I*(XHI-XLO ) 

DRVX=0. 

DO  30  J=1 , NDEGM1 

DRVX  =(DRVX+DC0EFF(NDEG-J+1 ))*X 
30  CONTINUE 

DRVX=DRVX*DCOEFF( 1  ) 

TEST=DHVXLO*DRVX 

IF(TEST.GT.O)  GO  TO  40 
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40 


IERR0R-1 

RETURN 

CONTINUE 

IERRORK) 

RETURN 

END 
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Q  ********************************************************************** 


c  *  * 

C  *  SUBROUTINE  ' PLOT1 '  -  CALLING  PROGRAM  NAME:  'ASSAY'  * 
C  *  * 
C  *  WRITTEN  MAY  1981  -  J.J.KNUDSEN,  LAIR,  DIV.OF  RESEARCH  SUPPORT,  * 
C  *  ANALYTICAL  CHEMISTRY  GROUP.  * 
C  *  * 
C  *  LANGUAGE  -  DATA  GENERAL  FORTRAN  V  * 
C  *  * 
C  #  PURPOSE  OF  SUBROUTINE  * 
C  * .  * 


C  *  SET  UP  A  TWO  DIMENSIONAL  ARRAY  CALLED  'GRAPH'.  THE  IMAGE  GRID  HAS  * 
C  *  'H'  ROWS  AND  'W  COLUMNS.  THE  ARRAY  IS  INITIALLY  FILLED  WITH  * 
C  *  BLANKS.  THE  BOUNDARIES  OF  THE  ARRAY  ARE  THEN  DEFINED  WITH  APPP.OP-  * 
C  *  RIATE  SYMBOLS  AND  TICK  MARKS  AND/OR  GRID  LINES  ADDED.  GRID  LINES  * 
C  *  ARE  OPTIONAL.  * 


C  *  * 

C  *  SUBROUTINE  FORMAL  PARAMETERS  * 

C  * . * 

C  *  INTEGER  * 

C  * -  * 

C  *  GRAPH  -  TWO  DIMENSIONAL  ARRAY  SETUP  WITH  GRID  LINES  IN  THIS  SUB-  * 

C  *  ROUTINE.  * 

C  *  * 

C  *  SUBROUTINE  SYMBOL  DEFINITIONS  (EXCLUDING  FORMAL  PARAMETERS)  * 

C  * .  * 


C  *  INTEGER  VARIABLES 


C  * -  * 

C  *  W  -  VARIABLE  SET  =  TO  WIDTH  OF  GRAPH  (100  SPACES)  * 
C  *  H  -  VARIABLE  SET  =  TO  HEIGHT  OF  GRAPH  (90  SPACES)  * 
C  *  BLANK  -  VARIABLE  SET  -  TO  ALPHAMERIC  CHARACTER  '  '  * 
C  *  PLUS  -  VARIABLE  SET  «  TO  ALPHAMERIC  CHARACTER  '♦'  * 
C  *  MINUS  -  VARIABLE  SET  -  TO  ALPHAMERIC  CHARACTER  * 
C  *  DOT  -  VARIABLE  SET  *  TO  ALPHAMERIC  CHARACTER  ’.'  * 
C  *  I  -  VARIABLE  USED  AS  A  COUNTER  * 
C  *  J  -  VARIABLE  USED  AS  A  COUNTER  * 


* 


SPECIAL  INSTRUCTIONS 


* 
* 

* _ _  * 

*  USE  /X  SWITCH  AT  COMPILATION  TIME  TO  DELETE  GRID  LINES  FROM  GRAPH  * 

*  ARRAY.  * 

*  * 
**#************»•#***#***«»*****»•»«****•**•»«•*»******•*•****•*»*•»•* 


SUBROUTINE  PL0T1  (GRAPH) 

INTEGER  W,H, GRAPH (0:90,0: 100) .BLANK, PLUS, MINUS, DOT, i ,J 

BLANK’’  ' 

PLUS’* ♦' 
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MINUU= ’ - ‘ 
!)0T- '  .  ' 


W=1'J0 


J  'II 

100  1*0 

200  GUAPH( J, I ) “BLANK 
1=1  +  1 

IP(I.LE.W)  GO  TO  200 
J=J-1 

IF(J.GE.O)  GO  TO  100 

CONTINUE 

J=H 

250  K=0 

300  1=0 

400  GRAPH (J,I+K)=PLUS 

IF(K.EQ.W)  GO  TO  600 
500  1-1+1 

X  I F  ( J .  EQ .  50 .  OR .  J .  EQ .  0 )  GO  TO  530 

X  GO  TO  540 

530  GRAPH(J,I+K)=MINUS 
540  1F(I.LT.9)  GO  TO  500 

K=K+10 
GO  TO  300 

600  IF(J.EQ.O)  GO  TO  1000 
L=  1 

620  J-J-1 

K=0 

640  1=0 

645  CONTINUE 

X  IF(l*K.EQ.O.OR. I*K . EQ. 100)  GO  TO  650 
X  GO  TO  660 

650  GRAPH (j , J*K )-DOT 

660  IF(I.EQ.IO)  GO  TO  800 
I-I+1 
K-1 0 

GO  TO  645 
800  L=L+1 

IF(L.LT.5)  GO  TO  620 
J=J-1 

IF(J.GE.O)  GO  TO  250 

1000  RETURN 

END 
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ooooooooooooooooooooooo 


0  »**«***•***»***•»*»*******•*•**•*•»»***##*#****«•**•*»•*»*•»*•**»*•»»* 


c  *  * 

C  *  SUBROUTINE  'PLOT 2'  -  CALLING  PROGRAM  NAME:  'ASSAY'  * 

C  *  * 

C  *  WRITTEN  MAY  1981  -  J.J.KNUDSEN,  LAIR,  DIV.OF  RESEARCH  SUPPORT,  * 

C  *  ANALYTICAL  CHEMISTRY  GROUP.  * 

C  *  * 

C  *  LANGUAGE  -  DATA  GENERAL  FORTRAN  V  * 

C  *  * 

C  *  PURPOSE  OF  SUBROUTINE  * 

*  _______  _  __________  * 


*  THIS  SUBROUTINE  PUTS  'SYMBOL'  INTO  PLOTTING  ARRAY  'GRAPH'  CREATED  * 

*  BY  SUBROUTINE  ' PL0T2 '  AT  A  POINT  CORRESPONDING  TO  (X,Y)  IN  THE  * 

*  IMAGE  PLANE  FOR  BOTH  INPUT  RAW  DATA  FOR  STANDARD  CURVE  AND  GEN-  * 

*  ERATED  REGRESSION  EQUATION.  THE  RANGE  OF  COORDINATES  IS  ASSUMED  * 

*  TO  BE  (XMIN  TO  XMAX,  YMIN  TO  YMAX)  AS  DEFINED  BELOW.  THE  PAIR  OF  * 

*  REAL  COORDINATES  (X,Y)  IN  THE  IMAGE  PLANE  ARE  CONVERTED  TO  A  PAIR  * 

*  OF  INTEGER  SUBSCRIPTS  (l,J)  USING  THE  DISCRETIZATION  FORMULA:  * 


*  * 

*  I  *  (X  -  XMIN/XMAX  -  XMIN)  *  N  * 

*  J  -  (Y  -  YMIN /YMAX  -  YMIN)  *  N  * 

*  # 

*  WHERE:  N  -  DESIRED  NUMBER  OF  DISCRETIZATION  INTERVALS  * 

*  • 


*  SUBROUTINE  FORMAL  PARAMETERS  * 

*  _ _ _ _ _  * 


*  INTEGER  * 

*  _ _  * 


*  GRAPH  -  TWO  DIMENSIONAL  ARRAY  SET  UP  BY  'PLOT1 ' .  VALUES  OF  INPUT 

*  STANDARD  CURVE  DATA  POINTS  INSERTED  IN  THIS  ARRAY  AS  WELL 

*  AS  50  EQUALLY  SPACED  POINTS  OF  REGRESSION  EQUATION  GEN- 

*  ERATED  BY  'ASSAY'. 

*  N  -  NUMBER  OF  STANDARD  CURVE  DATA  POINTS  ENTERED  IN  MAIN  PROGRAM 
C  *  NDEG  -  DEGREE  OF  REGRESSION  EQUATION  GENERATED  BY  'ASSAY' 

C  *  A  -  VARIABLE  SET  *  TO  0  IF  'MAXABS'  IN  MAIN  PROGRAM  <-  1.0 
C  *  SET  *  TO  1  IF  'MAXABS'  IN  MAIN  PROGRAM  >  1.0 

C  *  REAL 
C  * - 

C  *  DAT  -  TWO  DIMENSIONAL  ARRAY  CONTAINING  STANDARD  CURVE  DATA  POINTS 

C  * 

C  *  DOUBLE  PRECISION 
C  * . 

C  *  COEFF  -  VECTOR  OF  LENGTH  6  CONTAINING  REGRESSION  COEFFICIENTS 

C  * 

C  *  SUBROUTINE  SYMBOL  DEFINITIONS  (EXCLUDING  FORMAL  PARAMETERS) 

C  * . 

C  *  INTEGER  VARIABLES 
C  * . 

C  *  SYMBOL  -  VARIABLE  CONTAINING  SYMBOL  TO  BE  USED  IN  PLOTTING  INPUT 


* 

« 

* 


* 

* 

* 


* 

« 

» 

* 


* 
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c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 

c 

r 

c 

c 

c 

c 

c 

c 

c 

c 


*  STD  CURVE  DATA  POINTS  AND  REGRESSION  EQUATION  POINTS  * 

*  SET  -  TO  '0*  FOR  STD  CURVE  DATA  POINTS  * 

*  SET  «T0  ’(a)'  FOR  REGRESSION  EQUATION  DATA  POINTS  * 

*  W  -  WIDTH  OF  ARRAY  'GRAPH '  * 

*  H  -  HEIGTH  OF  ARRAY  'GRAPH'  * 

*  I  -  VARIABLE  TO  CONTAIN  'X-AXIS*  INFORMATION  FOR  PLOTTING  * 

*  J  -  VARIABLE  TO  CONTAIN  'Y-AXIS*  INFORMATION  FOR  PLOTTING  * 

*  NDEGP1  -  DEGREE  OF  REGRESSION  EQUATION  ♦  1  (LOOP  CONTROL)  * 

*  * 


*  DOUBLE  PRECISION  VARIABLES  * 

*  _ _  * 

*  XLO  -  LOWER  LIMIT  OF  'STANDARD  VOLUME  EQUIVALENT'  * 

*  XHI  -  UPPER  LIMIT  OF  'STANDARD  VOLUME  EQUIVALENT'  * 

*  X  -  VARIABLE  USED  IN  SOLUTION  OF  REGRESSION  EQUATION  FOR  PLOTTING  * 

*  PURPOSES.  'X-AXIS'  INFORMATION.  * 

*  Y  -  VARIABLE  USED  IN  SOULTION  OF  REGRESSION  EQUATION  FOR  PLOTTING  * 

*  PURPOSES.  'Y-AXIS'  INFORMATION.  * 

*  XMIN  -  MINIMUM  VALUE  OF  'X'.  * 

*  XMAX  -  MAXIMUM  VALUE  OF  'X'.  * 

*  YMIN  -  MINIMUM  VALUE  OF  'Y'.  * 

*  YMAX  -  MAXIMUM  VALUE  OF  'Y'.  SET  -  TO  1.0  IF  A  -  0,  SET  =  TO  * 

*  2.0  IF  A  -  1.  * 

*  * 


********************************************************************** 


SUBROUTINE  PLOT 2  (GRAPH,N,DAT,NDEG,COEFF,A) 

INTEGER  GRAPH(0:50, 0:100), SYMBOL, I, J,W,H,NDEG,NDEGP1 ,N,A 
REAL  DAT (30,3 ) 


DOUBLE  PRECISION  XHI , XLO, C0EFF(6 ) ,X, Y, XMIN , XMAX, YMIN , YMAX 


IF(A.EQ.O)  GO  TO  20 
YMAX=2.0 
GO  TO  30 
20  YMAX- 1.0 

30  XMAX -1.0 

YMIN-0.0 
XMIN -0.0 
XLO-O.O 
XHI-1 .0 


H-51 

W-100 

SYMBOL- 'O' 

DO  100  K-1 ,N 

I-INT(DAT(K,3)*(DBLE(FL0AT(H)))) 

J«INT(((DAT(K,2)-YMIN)/(YMAX-YMIN))*(DBLE(FL0AT(W)))) 
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1 

1 

1 


I »50-I 


IF((I.GE.0)  .AND. 

(I.LE.H)  .AND. 

(J.GE.O)  .AND. 

(J.LE.W)  )  GO  TO  75 

GO  TO  100 
75  GRAPH ( I , J ) “SYMBOL 

100  CONTINUE 

SYMBOL-’ O' 

NDEGP1 “NDEG+1 
DO  200  K-0,50 

X*XL0+O.O2#K*(XHI-XL0 ) 

Y-0. 

DO  150  L=1 ,NDEG 

Y-(Y+C0EFF(NDEGP1-L+1  ))*X 
150  CONTINUE 

Y=Y+C0EFF(1 ) 

I-INT(X*(DBLE(FLOAT(H)))) 

j*int(( (y-ymin)/ (ymax-ymin))*(dble(float(w) )) ) 
I-50-I 

IF((l.GE.O)  .AND. 

1  (I.LE.H)  .AND. 

1  (J.GE.O)  .AND. 

1  (J.LE.W)  )  GO  TO  175 

GO  TO  200 

175  GRAPH ( I , J ) “SYMBOL 


200  CONTINUE 

RETURN 
END 
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